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Ion acceleration driven by superintense laser pulses is attracting an impressive and 
steadily increasing effort. Motivations can be found in the potential for a number of 
foreseen applications and in the perspective to investigate novel regimes as far as avail- 
able laser intensities will be increasing. Experiments have demonstrated in a wide range 
of laser and target parameters the generation of multi-MeV proton and ion beams with 
unique properties such as ultrashort duration, high brilliance and low emittance. In this 
paper we give an overview of the state-of-the art of ion acceleration by laser pulses as 
well as an outlook on its future development and perspectives. We describe the main 
features observed in the experiments, the observed scaling with laser and plasma pa- 
rameters and the main models used both to interpret experimental data and to suggest 
new research directions. 
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I. INTRODUCTION 



More than half a century ago, Veksler (1957) intro- 
duced the concept of "coherent acceleration" of particles 
as a mechanism in which the accelerating field on each 
particle is proportional to the number of particles being 
accelerated, in contrast to traditional techniques. Addi- 
tional elements in Veksler's vision of future accelerators 
included the automatic synchrony between the particles 
and the accelerating field, the localization of the latter 
inside the region where the particles are, and the pro- 
duction of quasi-neutral groups with large numbers of 
energetic particles. 

These features are realized in the acceleration of ions 
from plasmas produced by intense laser pulses. There, 
as a very general description, strong electric fields are 
generated by a collective displacement of a large number 
of electrons, and such electric fields accelerate ions until 
charge neutrality is restored and ions move together with 
electrons in a ballistic way. 

Before the year 2000, ions having energies up to several 
MeVs had been observed in several high-intensity laser- 
matter interaction experiments and for different targets, 



including thick solid targets (Beg et al. 1997 Fews et al. 



1994 



Gitomer et al.\ 1986 and references therein), gas 



jets ( 
refere 


Krushelnick et al. 


1999 


Sarkisov et al. 


1999 


and 


mces therein) and sub-micrometric clusters ( 


Ditmire 



et al. 1997 1999 and references therein). Common to 
these experiments was the rather isotropic ion emission 
and the resulting low brilliance, making these configura- 
tions not attractive as ion accelerators for applications. 



In 2000 three experiments (Clark et al. 2000a Mak 



simchuk et al. 2000 Suavely et al. 20001 indepen- 



dently reported the observation of an intense emission 
of multi-MeV protons from solid targets, either metal- 
lic or plastic (CH), of several microns thickness irradi- 
ated by high intensity laser pulses. The basic set-up 
of these experiments in shown in Fig{T] The laser in- 
tensities, number of protons and maximum ion energy 



observed for the three experiments were, respectively. 
3x10 ^^ W cm-2 
2000|), 5 X 10 ^^ 

1020 



et al. 



> 10^ and 1.5 MeV ( [Maksimchuk al. 
W cm-^ - 10^2 and 18 MeV (IClark 



2000a), and 3 



58 MeV (Snavely et al. 



W cm-^ - 2 X lO^^^Mid 

2000). Fig 2^ shows the spectrum 



2|show 

" pooo 



of ions observed by Snavely et al. (12000). The protons 
were emitted as a quite coUimated beam in the forward 
direction with respect to the laser pulse propagation and 
were detected at the rear side of the target, opposite to 
the laser-irradiated surface. 

The emission of protons from metallic targets whose 
chemical composition does not include hydrogen may 
sound surprising, but it was already clear from previous 
experiments that protons originated from impurities, i.e. 
thin layers of water or hydrocarbons which are ordinarily 
present on solid surfaces. In experiments performed both 
with "long" , nanosecond pulses ( Gitomer et al.\ 1986', and 
references therein) and "short" (sub-) picosecond, high 
intensity pulses 



Fews et al.\ |1994 



Beg et al.\ p97| [Clark et al.\ |2000b 



protons and heavier ions were com- 
monly detected in the backward direction (i.e. towards 
the laser) with a broad angular distribution and their 
origin was interpreted in terms of acceleration in the ex- 
pansion of the hot laser-produced plasma at the front 
(laser-irradiated) side of the target. The characteristics 
of the forward proton emission in the new experiments, 
such as the high degree of coUimation and laminarity of 
the beam, were much more impressive. 

These findings generated an enormous interest both in 
fundamental research and in the possible applications. 
For these latter, the most relevant and peculiar feature 
of multi-MeV ions is the profile of energy deposition in 
dense matter. Differently from electrons and X-rays, pro- 
tons and light ions deliver most of their energy at the end 
of th eir path (Fig|3l), at the so-called Bragg peak (Knollj 



2010 Ziegler et al. 2008 ). The physical reason is that the 
energy loss is dominated by Coulomb collisions for which 
the cross section strongly grows with decreasing energy, 
so that the stopping process becomes progressively more 




FIG. 1 (Color online) Artist's view of a typical experiment 
on proton emission from laser-irradiated solid targets. 
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FIG. 2 Proton energy spectrum from tlie rear side of a 100 fim 
solid target irradiated by a 423 J, 0.5 ps pulse at normal inci- 
dence, corresponding to an intensity of 3 x 10^" W cm~^. The 
integrated energy of protons indicates a conversion efficiency 
of ~ 10% for protons above 10 MeV. Reprinted figure with 
permission from |Snavely et aZ.| ( [2000[ ), Phys. Rev. Lett. 85, 
2945. Copyright 2013 by the American Physical Society. 



and more efficient. This property makes protons and ions 
very suitable for highly localized energy deposition. The 
applications that were proposed immediately after the 
discovery of multi-MeV proton acceleration included ion 
beam cancer therapy, laser triggering and control of nu- 
clear reactions, production and probing of warm dense 
matter, "fast ignition" of Inertial Confinement Fusion 
targets and injectors for ion accelerators. These foreseen 
applications are reviewed in Sec|Vj As a particularly in- 
novative and successful application, ultrafast probing of 
plasmas by laser-driven proton beams will be described 
in SecfOI 

While the potential for applications was apparent, the 
details of the physics behind proton acceleration were not 
clear. A debate arose on the actual location of the region 
where the protons were accelerated and, consistently, on 



the mechanism driving acceleration. Clark et al. (2000a 



and Maksimchuk et al. (2000) suggested that protons 



were accelerated at the front side of the target, cross- 
ing the latter and being detected on the opposite side. 



In contrast, Snavely et al. (2000) provided evidence that 



protons were accelerated at the rear side (see also Hatch 
ett et al. ( 2000[ )). To support the interpretation of these 
latter experiments (performed at the Petawatt facility of 
Lawrence Livermore National Laboratory, USA) the so- 
called Target Normal Sheath Acceleration (TNSA) model 
was introduced bylWilks aL|(|2001|). Briefiy, TNSA is 



driven by the space-charge field generated at the rear 
surface of the target by highly energetic electrons accel- 
erated at the front surface, crossing the target bulk, and 
attempting to escape in vacuum from the rear side. The 
basic theory of TNSA and related models are described 



in detail in Sec. Ill Most of the experiments later per- 
formed on proton acceleration by laser interaction with 
solid targets have been interpreted in terms of the TNSA 



framework (Sees. III. A and III.B I that has also guided 
developments towards source optimization by target en- 



gineering (Sec. III.E). 

A major request for several of the foreseen applications 
is an increase of the energy per nucleon up to the hun- 
dred of MeVs and beyond. The next generation of laser 
facilities should allow intens ities higher than the cu rrent 



record of - 10^^ W cm"^ (Yanovsky et al. 2008), but 



at present it is neither guaranteed that the ion energy 
scaling observed so far will be maintained at such ex- 
treme intensities nor that TNSA will be still effective. An 
analysis of proton acceleration experiments performed up 
to 2006 suggests a ^ (/A^)^/^ scaling of proton energy 
with laser intensity / and w avelength A up to values of 
/A2 = 3 X 10^0 W cm" ^ 



"^/im^ 



(Borghesi et al. 



2006 



Fuchs 



eFaH|2006b| . FigjI] summarizes such data, together with 



more recent results obtained with Ti:Sa-based, ultrashort 
(tens fs) pulses, exhibiting a ^ IX^ scaling. Measure- 



ments by Robson et al. (20071 at energies up to 400 J, 



pulse durations between 1 and 8 ps and intensities up to 
6 X 10^" W cm~'^ suggested a lower scaling. It is also of 
crucial importance to establish the most relevant scaling 
parameters as well as improve or optimize beam emit- 
tance, brilliance and monoenergeticity for specific appli- 
cations. For instance TNSA-generated proton beams are 



highly laminar and have very low emittance (Sec. III.B ) 



but the energy spectrum is ordinarily broad and thus not 
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FIG. 3 Example of the profile of energy deposition of protons 
and C ions in water, compared to those of electrons, X- and 
7-rays, and neutrons. Protons and C ions profiles are charac- 
terized by the Bragg peak at the end of the path. The quan- 
tity plotted is the relative dose, i. e. the energy absorbed pe r 
unit mass. Reprinted figure from Amaldi and Kraft (20051, 
Rep. Prog. Phys. 68, 1861. By permission from Institute of 
Physics Publishing (2013). 
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FIG. 4 (Color online) Maximum proton energy from laser- 
irradiated solid targets as a function of the laser irradiance 
and for three ranges of pulse durations, and with additional 
data (labels "Lund" and "Saclay") reporting later experi- 
ments up to 2008. Two trendlines are overlayed, the shal- 
lower one corresponding to a ~ I^^^ dependence, and the 
steep er one to a sca l ing pr oportional to I. Reprinted figure 
from Borghesi et al. (20081, Plasma Phys. Contr. Fusion 50, 



124040. By permission from Institute of Physics Publishing 
(2013). 



optimal for most applications. 

These issues motivate the search for other ion accelera- 
tion mechanisms. The latter include concepts which have 
been explored previously in different contexts (e.g. astro- 
physics), such as Radiation Pressure Acceleration (RPA) 
and CoUisionless Shock Acceleration (CSA). Other pro- 
posed schemes exploit the potential both of advanced 
target engineering and of nonlinear "relativistic" optical 
effects in plasmas, such as ion acceleration in ultrathin 
solid targets which become transparent to intense laser 
pulses (Break-Out Afterburner, BOA) or involving low- 
density targets. The basic physics of these mechanisms 
and the related experimental work, still in a preliminary 



stage with respect to TNSA, will be described in Sec IV 



The development of advanced acceleration schemes is sus- 
tained by the continuous trend towards laser pulses of 
higher intensity and larger energy. A detailed account of 
the many running projects and developing facilities re- 
lated to optical and infrared lasers (A ~ 1 ^m) is given 
in Sec. II of Di Piazza et al. (2011 ). Progress towards CO2 
lasers (A ~ 1 /xm) having Terawatt power (Haberberger 



et al.\ 2010[ ) is also of growing interest for ion accelera- 
tion, as discussed in Secs lIV.ATT] and [TV.B| 

On the theoretical side, the interpretation of experi- 



ments has revitalized classic and often controversial prob- 
lems of plasma physics such as plasma expansion into 
vacuum and the formation of coUisionless sheaths, at the 
basis of most of TNSA models, as well as other general 
physics models such as the motion of relativistic moving 
mirrors, a concept already discussed in the original work 



on special relativity by Einstein ( 1905 ) , which serves as 



a model for RPA. Simple, analytically affordable models 
are extremely useful to understand the basic accelera- 
tion mechanisms and, in particular, to provide scaling 
laws which may give directions for further developments. 
Reference models also highlight the several connections 
with other fields, such as the physics of discharges, of 
ultracold plasmas, and of particle acceleration in astro- 
physics. The theoretical discussions in Secs |III| and |IV] 
are, to a significant extent, based on this approach. 

Beyond simple modeling a rich and complex dynamics 
of laser- plasma interaction and ion acceleration, involving 
collective and self-organization effects, is apparent. Un- 
folding such dynamics requires the use of self-consistent 
electromagnetic, kinetic simulations. To this aim, the 



particle-in-cell (PIC) method (Sec II. D ) is by far the 
mostly used approach. Large-scale, multi-dimensional 
PIC simulations running on parallel supercomputers are 
an effective support for the design and interpretation 
of laser-plasma acceleration experiments, although fully 
"realistic" simulations in three spatial dimensions and 
for actual laser and target parameters are most of the 
time still beyond computational possibilities. These lim- 
itations further motivate the development of complemen- 
tary, reduced simulation models. These issues are further 
discussed in Secs lII.DI and lIII.C.4l 

A comprehensive description of laser-plasma dynam- 
ics if far beyond the scope of the present work and 
may be found in recent books and reviews dGibboKl 



2005b Mourou et al. 



2006 Mulser and Bauer 20101. 



In Sees II. A and II. B we only describe a few basic is- 
sues of relevance for the understanding of ion accelera- 
tion mechanisms. These latter are first introduced in a 
compact way in Sec|II.C[ leaving a detailed discussion in 



the following Secs|III|and IV 



M. LASER ION ACCELERATION IN A NUTSHELL 

A. Laser interaction with overdense matter 

In the present work we mostly refer to ion acceleration 
occurring in the interaction with solid targets, for which 
the electron density Ue greatly exceeds the so-called crit- 
ical or cut-off density. 



= 1.1 X 10 



21 



A 



1 /im 



(1) 



The condition Ue = ric is equivalent to ujp — u) where 
LOp = {'innee'^ /nieY^'^ and w = 2t:c/\ are the plasma 
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and laser frequencies, respectively. Since the linear re- 
fractive index of the plasma is n = (1 — Wp/w^)^/^ = 
(1 — rie/nc)^/^, in the rig > ric "overdense" region n has 
imaginary values and the laser pulse cannot propagate. 
All the laser-plasma interaction occurs either in the "un- 
derdense" region rig < ric or near the "critical" surface 
at which Ue — Uc- 

Relativistic effects make the refractive index nonlinear. 
Qualitatively speaking, the "relativistic" refractive index 
describing the propagation of a plane wave whose vector 
potential is A = A(x,t) is obtained from the linear ex- 
pression by replacing the electron mass with the quantity 
TOe7, where the relativistic factor 7 is given by 



7 



(2) 



a = eA-lmeC^ , and the angular brackets denote an aver- 
age over the oscillation period. The parameter ao is the 
commonly used "dimensionless" amplitude related to the 
laser intensity / bj|^ 



ao = 0.85 



1/2 



1018 W cm" 



(3) 



where we used / = cl^E"^^ /Air to relate the electric field 
E = —{l/c)dA/dt to the laser intensity /. 

The nonlinear, "relativistic" index p^l = 
(1 — 7ie/(7rtc))^^^ becomes imaginary when rig > "/ric, 
showing an increase of the cut-off density for a plane 
wave: this effect is known as relativistic self-induced 
transparency or, briefly, relativistic transparency. How- 
ever, the problem of laser penetration inside a plasma 



is not such trivial ( 


Cattani et al. 


2000 


Goloviznin and 


Schep 


2000 


Shen and Xu 


2001 


) because of both the 



nonlinearity in the wave equation and the self-consistent 
modification of the plasma density profile due to radi- 
ation pressure effects. These latter may be described 
via the ponderomotive force (PF)j^ In an oscillating, 
quasi-monochromatic electromagnetic field described by 
a dimensionless vector potential a(r, t) whose envelope 
is sufficiently smooth in space and time, the relativistic 
PF is [see e.g. Bauer et al. { 1995| ) and Mulser and Bauer 
([2OIOI)] 



-mcC^V{l + (a) 



xl/2 



^ Consistently with our definitions, given the value for I the peak 
value of the dimensionless vector potential of the plane wave will 
be given by ao for linear polarization and by ao/V2 for circular 
polarization. 

^ Throughout the present review we refer to the ponderomotive 
force as the slowly- varying, effective force describing the cycle- 
averaged motion of the "oscillation center" of a charged particle 
in an oscillating non-uniform field, over a time scale longer than 
the oscillation period. "Fast" oscillating components are not 
included in the definition of ponderomotive force here adopted. 



For a plane wave impinging on an overdense plasma, the 
resulting PF, more effective on the lightest particles, is 
in the inward direction and tends to push and pile up 
electrons inside the plasma. Ponderomotive effects will 
be further discussed below (see Sees II.B|II.C ) . 

In a multi-dimensional geometry, a laser pulse of fi- 
nite width may produce a density depression around the 
propagation axis also because of ponderomotive pushing 
of the electrons in the radial direction. Jointly with the 
relativistic effect and target expansion driven by elec- 
tron heating, this mechanism may lead to a transition to 
transparency as soon as the electron density drops the 
cut-off value (Fuchs et al. 1999). Investigations of ion 



acceleration in the transparency regime are described in 
Sec lRTCl 

The penetration of the laser pulse depends not only 
on the electron density but also on the target size when 
the latter becomes close to, or smaller than one wave- 
length. As a simple but useful example, the nonlinear 
transmission and reflection coefficients can be calculated 
analytically for a sub-wavelength foil modeled as a Dirac 



delta- like density profile ( Macchi et al. 2009c Vshivkov 



et al. 1998), obtaining a transparency threshold 



Ur, A 



(5) 



where £ is the thickness of the foil. This formula has some 
interest for the interaction with ultrathin foil targets (see 
Secs |IV.A.2pv:C ). 



B. Hot electrons 

Since the laser pulse cannot penetrate into solid density 
regions, the absorbed energy is there transported mostly 
by energetic (commonly named either "hot" or "fast") 
electrons which may be generated during the interaction 
by several mechanisms. By "hot" electrons in the present 
context one typically refers to relativistic electrons whose 
energy is of the order of the cycle-averaged oscillation 
energy in the electric field of the laser in vacuum. 



£p = rrieC^i-y - 1) 



\/l + a§/2 



1 



(6) 



(4) where Eq. ^ has been used. Expression ([6|) is also called 



the "ponderomotive" energy (Wilks et al. 119921). Hot 
electrons penetrating into solid targets have been ob- 
served and characterized in several experiments at very 
high intensities and for different interaction condition^ 
and play a fundamental role in applications such as laser- 
driven photonuclear physics and fast ignition of fusion 



See e g . |Chen et al.\ ||2009a|l; [Key et al.\ (|l998t; [Nilson et al.\ 
(2010'); 'Tanim oto et a;.| | |2009l l; and [Wharton et a;.| | |1998| l arid 
references therein. 
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targets. Moreover, as it will be discussed in Sec. |III[ in 
most of the experiments reported so far, acceleration of 
protons and heavier ions is driven by hot electrons. 

The process of hot electron generation turns out to be 
complex and, possibly, not completely understood yet. 
A complete account of past and ongoing research on 



cleaning techniques may allow to minimize prepulse ef- 



the topic may be found in recent books ( Gibbon 2005b 



Mulser and Bauer, 2010) and in a vast experimental and 



theoretical literature. Here we give a very basic discus- 
sion at a qualitative level, focusing on those aspects which 
are most essential and relevant to ion acceleration. 

At the surface of an overdense plasma, electrons are 
driven by the Lorentz force — — e(E + v x B/c) which 
include both the incident and reflected laser pulses and 
the self-generated fields. As a necessary condition for 
the efficient generation of hot electrons near the critical 
surface must have an oscillating component directed 
along the density gradient Vrig. This is the case for 
the well-known resonance absorption where the condi- 
tion E • Vrie ^ is necessary to drive resonant plasma 
oscillations which in turn accelerate electrons. In a plane 
geometry such condition requires oblique incidence and 
P-polarization of the laser pulse, and that absorption is 
sensitive to the density scalelength L„ = n^/l'Vu^l be- 
cause the driving force is evanescent in the resonance 
region. 

In a sharp boundary plasma where L„ <C A, absorp- 
tion and heating may arise because electron motion is 
not adiabatic, as electrons are driven from the region of 
strong fields to the evanescence region in a time shorter 
than 2n /(jj, so that the cycle average — e (E • v) may not 
cancel out. Thus, short duration and high intensity laser 
pulses favor electron heating because the hydrodynamic 
expansion has both not sufficient time to wash out sharp 
density gradients ot it is dominated by the strong pon- 
deromotive force that steepens the density profile. 

At this point it is worth to remind that in most high- 
intensity experiments the main interaction pulse is pre- 
ceded by prepulse^ which cause early plasma formation 
and expansion, so that the short pulse interaction does 
not occur with a sharp-boundary, solid-density plasma. 
However, profile steepening at the critical surface will 
be still effective, thus one may expect the interactions 
still occurs with a sharp density profile, having with a 
lower density jump with respect to a solid target. Oc- 
casionally, "preplasma" formation may also allow addi- 
tional electron acceleration mechanisms to take place in 
the underdense plasma region (Esarey et al. 2009), pos- 



sibly leading to electron energies much higher than given 
by Eq. In more recent experiments, advanced pulse 



fects (see Sec III.E). 



1. Heating models 

A popular electrostatic model of electron heating at 
a step-boundary plasma has been proposed by |Brunel| 
(1987). In this model, electrons in this system are 



dragged out of the surface of a perfect conductor by an os- 
cillating "capacitor field" , extending on the vacuum side, 
and representing the P-component of the incident plus 
refiected laser electric field. Electrons are considered to 
be "absorbed" when, after having performed about half 
of an oscillation on the vacuum side, they re-enter the 
target there delivering their energy, which is of the or- 
der of the oscillation energy in the external field|^ The 
model thus accounts in a simplified way for the pulsed 
generation (once per cycle) of hot electrons directed into 
the target and having an energy, roughly speaking, close 
to the "vacuum" value This simple model is not 

self-consistent because, for instance, the capacitor field is 
assumed to vanish inside the target, implying the pres- 
ence of a surface charge density. Nevertheless, following 



Mulser et al. (2001) it is possible to provide a "minimal" 



ID model, still in the capacitor approximation, where the 
electrostatic field is calculated self-consistently and an ac- 
celeration of electron bunches similar to that inferred by 
Brunei is apparent. We take the electric field as the sum 
of the electrostatic and a driver fields, e.g. = Ep + Ed 
where Ed = Ed{t) sinujQt with Edit) a suitable temporal 
envelope, a step-boundary density profile = noQ{x) 
{Z = 1 for simplicity) and a "cold" plasma, i.e. we ne- 
glect thermal pressure. Following these assumptions we 
write Maxwell's equations for the electrostatic field and 
Euler's equation for the electron fluid having velocity 



dxEe = 4:TTp = 47re[noO(a;) 



dvx 
~dt 



= {dt + Vxdx)vx = 



-iE,+Ed). 



(7) 
(8) 

(9) 



Switching to Lagrangian variables xq and ^ = £,{xo,t) 
defined by x = xq + d£,/dt = Vx, a. straightforward 
calculation along with the constraint of Eg being contin- 
uous at a; = a;o -|- C = yields the following equations 
of motion describing electrostatic, forced oscillations of 
electrons across a step-like interface: 



dt^ 



-UJ^Xo ~ eEd/lTLe 



(xo + e > 0) 
[xo + e < 0) 



(10) 



* In general the main pulse is preceded by both short pulses of ^ This effect is also commonly refereed to as "vacuum heating", 
similar duration as the main pulse, and by a much longer pedestal See [Gibbon] | |2005b| l for a discussion also on the origin of the 

due to amplified spontaneous emission. name. 
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FIG. 5 Numerical solution of the electrostatic "plasma sheet" 
model based on Eqs.( |10[ ) plus the exchange of initial position 
for crossing plasma sheets (see text for details). The trajec- 
tories of a limited number of sheets (one over twenty) in the 
{x, t) plane are shown. The driver field has the profile of an 
evanescent wave with peak amplitude Q.brrieUJc/ e in vacuum 
and a sin^(7ri;/2r) rising front with r = 5T where T = 2-r /uj. 
A density ne/ric = (ujp/u!)^ = 5 is assumed. 



From Eqs.(lO) we see that electrons crossing the bound- 
xo + ^ < 0) feel a secular forc e uj'^xq leading 



ary {x 

to dephasing from Ej, and acceleration (iMulser et al. 



2001). Eqs.(10l can be solved numerically for a discrete 



but large ensemble of electron "sheets" (corresponding to 
a set of values of xo > 0), with the prescription to ex- 
change the values of xq for two crossing sheets to avoid 
the onset of singularity in the equations]^ Representative 
trajectories of electrons moving across the interface are 
found as in FigjSj Electrons whose trajectory extends in 
vacuum for half or one period of the driving field and 
then re-enter at high velocity inside the plasma are ob- 
served. Similar trajectories are found in electromagnetic 



and self-consistent simulations (Sec |II.BT2 ). 

Notice that the "cold" plasma assumption is consistent 
with the requirement that the external field should be 
strong enough to overcome the potential barrier which, 
in an equilibrium state, confines warm electrons inside 
a bounded plasma (such barrier corresponds to a Debye 



sheath, see also Sees II.C.l and III.C.l ). For cUp ^ uj and 
nearly total reflection, the laser field component normal 
to the surface has an amplitude E± ~ 2EQsm9, with 
Eq = {inl/cY^^ the amplitude in vacuum. The sheath 
field Es ~ Te/{e\D) = (47rnoTe)i/2^ ^hat the condi- 
tion Ej\_ > Es may be rearranged as 4(//c) sin'^ 9 > noTe- 
This implies (at non-grazing incidence) the radiation 
pressure to exceed the thermal pressure and thus to coun- 
teract the thermal expansion and to steepen the density 



profile, making the assumption of a step-like plasma more 
self-consistent. 

For 5'-polarization or normal incidence there is no com- 
ponent of the electric field perpendicular to the surface. 
However, for high intensities the magnetic force term 
becomes important and may drive electron oscillations 
along the density gradient also for normal incidence. This 



effect is commonly named as "J x B" heating ( Kruer and 



Estabrook 1985). By considering the driver capacitor 



field as a model for the magnetic force component, the 
related electron dynamics may still be described using 
the above outlined models, but with two significant dif- 
ferences. First, to lowest order the magnetic force oscil- 
lates at 2w, thus leading to the generation of hot electron 
bunches twice per laser period. Second, the oscillating 
component perpendicular to the surface vanishes for cir- 
cular polarization (and normal incidence), so that hot 
electron generation might be strongly suppressed in such 
conditions. In fact, the vector potential representing a 
plane, elliptically polarized field may be written as 



A(x,t) = 



_Mx) 

\/rT< 



z{y cosujt + ezsinujt), (11) 



Using B = V X A and pj^ = eA/c for the transverse 
momentum of electrons, the — e(v x B/c) force can be 
written as 



V 

e- X B 

c 



.^e^d^AHx) 



1 + 



1 



1 + e^ 



DS 2a;t )(12) 



showing that the oscillating component vanishes for cir- 
cular polarization (e = l)j^ 

The integral over x of Uefpx, where fp^ = (fx) = 
(— e(v X B/c)a;) is the steady ponderomotive force den- 
sity on electrons, equals the total radiation pressure on 
the target surface. For circular polarization and normal 
incidence we thus expect radiation pressure to push the 
target while electron heating is quenched. These condi- 
tions have been investigated in order to optimize radi- 
ation pressure acceleration of ions versus other mecha- 



nisms driven by hot electrons, see Sec IV. A 



2. Simulations, multi-dimensional effects and simple estimates 

A more quantitative description of laser absorption and 
hot electron generation requires numerical simulations. 
To address electromagnetic effects in his model |Brunel| 
( 1988 ) performed two-dimensional (2D) PIC simulations 
in a plane wave, oblique incidence geometry. Several later 
studies using ID simulations with the "boosted frame" 



technique (Bourdier 19831 are summarized and reviewed 



® This numerical implementation basically corresponds to the pi- 
oneering, elementary model of plasma simulation formulated by 
|Dawson| | |1962| . 



^ A more detailed analysis shows that electron heating is quenched 
when the parameter e exceeds some threshold value, see |Macchi| 
I et a/.| ( |2009b[ | and [Rykovanov et al.\\2008\ . 
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by Gibbon et al. (1999). The absorption degree of a P- 
polarized laser pulse is quite sensitive to the incidence 
angle and the density scalelength, with the latter vary- 



ing on the timescale of ion motion ( Gibbon 1994 1 yield- 



ing a time-dependent absorption. Experimental attempts 



(Batani et al. 


2010 Flacco et al. 


2008 


McKenna et al. 


2008 


) have been made to vary the density scalelength in 



order to increase absorption in hot electrons and conse- 



quently to enhance ion acceleration (see Sec III.El. Hot 
electron generation tends to become more efficient for 
lower plasma densities, and particularly close to the criti- 
cal density Uc, as it is observed that stronger coupling and 
volumetric heating occurs near the transmission thresh- 
old. A "near-critical" plasma may be either produced by 
the laser prepulse or by using a special target material, 



e.g. a low-density foam (see Sec IV.D). 

2D simulations reveal additional effects, as for instance 
the deformation of the plasma surface due to radiation 



pressure-driven "hole boring" (see also Sec IV. A. 1) that 



changes the local incidence angle ( Wilks et al.\ 119921), 
leading to increased absorption and providing a dynamic 
"funnel" effect collimating the electron flow inside the 
target (Ruhl et al. 1999). A similar dynamics occurs in 



microcone targets which have proved to be effective in 



enhancing hot electron generation (Gaillard et al. 2011 
[Nakamura eTd] |2009| [Sentoku et al.\ |2004| a^ refer 
ences therein). 

Absorption is also sensitive to small-scale surface de- 
formations, either self-generated or pre-imposed, so that 
the use of microstructures on the front target surface of 
has been also suggested as a way to enhance hot electron 
generation: see e.g. Klimo et al. (2011) and references 



therein. Another possible approach is the use of grating 
surfaces where the resonant excitation of surface plasma 



waves may also lead to very high absorption (Bigongiari 



et al. 2011 Raynaud aL 2007). 



The high sensitivity of hot electron generation to laser 
and plasma parameters partly accounts for data scatter 
and differences observed in the many experimental inves- 
tigations reported in the literature, with the above men- 
tioned prepulse effects bringing additional complexity. 
For these reasons, absorption values and characteristics 
of the hot electron distribution are often taken into ac- 
count in a phenomenological way. It has been often con- 
sidered acceptable to assume the hot electron distribution 
to be Maxwellian with a temperature Th given by Eq.([6| 
as a function of the laser irradiance. Fig[6]presents a col- 
lection of temperature measurements obtained for sub- 



picosecond pulses up to the year 2000 ( |Gibbon[ 2005b|; 
these data broadly support a scahng of Th as (/A^p^. 
The total fractional absorption in hot electrons 77/1 is usu- 
ally estimated to be in the 10% — 30% range, with ex- 
perimental indications of possibly quite higher values at 
ultra-relativistic intensities ( |Ping et al. 2008). An en- 
ergy flux balance condition such as rjhl — nuVhTh (with 
Vh — c a,t ultra-high intensities) may then be used to 
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FIG. 6 (Color online) Hot electron temperature as a function 
of irradiance from experiments of sub-ps laser-solid interac- 
tion. See table 5.2 in Gibbon ( 2005b[ ) for details on experi- 
mental parameters, diagnostic methods and references. The 
lines give scaling laws derived from different models [FKL: 



Forslund et al] ( (T977|), W: |Wilks et al.\ ( p92| , GB: [Gibbon 



and Bell 



(1992h, B: Brunei 



( 1987 i]T Reprinted figure from 
Gibbon ( 2005b|, Short Pulse Laser Interaction with Matter 
Imperial College Press). By permission from Imperial Col- 
lege Press and World Scientific (2013). 



estimate the "initial" density of hot electrons rih^ which 
usually is not larger than nc^ consistently with the ar- 
gument that Uh cannot exceed the density of the region 
where hot electrons are generated. 

Inside the target, the effective density may become dif- 
ferent from the above estimate for rih in particular condi- 
tions due to, e.g., the angular divergence of the electron 
flow or to electron refluxing effects depending on the elec- 
tron time of flight and recirculation and thus on the tar- 



get thickness (Mackinnon et al. 2002). Still one might 



roughly estimate the total number of hot electrons Nh 
by an energy balance relation Nh rihUh/Th where Ul 
is the energy of the laser pulse. The angular divergence 
6'div is also estimated from experiments to range between 



20 and 60 degrees and to increase with irradiance ( Green 



et al. 2008 and references therein), although such esti- 



mates might depend on the accuracy of the sheath field 



modeling (Ridgers et al. 2011). 



3. Hot electron transport in solid matter 

Transport of hot electrons in solid matter has been 
investigated extensively also because of its relevance for 



the electron-driven Fast Ignition scheme in ICF [see Free- 
man et al. ( 2006 ) for a survey] . Key issues characterizing 
this regime are the very high values of the currents and 
the effect of self-generated fields. From the above esti- 
mates it can be inferred that near the front surface of 
the target the current density jj^ = —erihVh associated 



9 



to hot electrons may reach values up to jh ~ ericC ~ 
4.8 X 10^^ A cm"^, corresponding to a total current of 
~ 15 MA over a circular spot of 10 /zm radius. This huge 
current must be locally neutralized by a return current 
jr such jh + jr — 0, otherwise either the electric field 
generated by the charge unbalance or the magnetic field 
generated by the free flowing current jh would be strong 



enough to stop the hot electrons ( Davies et at. 1997 Pas- 
soni et al. 2004 ) . The free, "cold" electrons contributing 



to the return current are either present as conduction 
electrons in metals or produced by field and coUisional 
ionization in insulators ( Tikhonchuk 2002 1 . Additional 



complexity is introduced by effects such as target heat- 
ing and hot electron refluxing, which have been inferred 



in several experiments (Bellei et al. 2010 Nilson et al. 



2011 Quinn et al. 2011 and references therein). Fil- 



amentation instabilities and dependence on the target 



material have also been extensively studied ( Fuchs et al. 



[20031 iMaSclossi aLH2006l [McKenna et aZ.l|20111 and 
references therein). Simulation models accounting for 
both coUisional effects and self-consistent generation of 
quasi-static fields are needed for quantitative investiga- 
tions!^ Finally, it is noticeable that at least a fraction 
of hot electrons propagate coherently through the target 
conserving the temporal periodicity of the driving force, 
i.e. as bunches with lo or 2w rate depending on incidence 
angle and polarization, as inferred by optical transition 



radiation measurements ( Popescu et al. 2005 ) 



C. Ion acceleration mechanisms 

In this Section we give an overview of ion accelera- 
tion mechanisms including both those proposed to ex- 
plain early experimental results in solid targets and those 
investigated later, either following inspiration from theo- 
retical work or testing novel target designs. Some of the 
mechanisms described below and the target regions where 
they are active are indicated in the cartoon of Figj?] Ion 
acceleration models will be described more in detail in 
Sees |III| and |I V| along with the most relevant experiments. 



1. Rear surface acceleration 



As outlined in Sec |II.B[ a very intense current of high- 
energy hot electrons may be generated at the front side 
of the target and eventually reach the rear side. There, 
as the hot electrons cross the rear side boundary and 
attempt to escape in vacuum at the rear side, the charge 
unbalance generates a sheath field Eg , normal to the rear 



See e.g. |Bell et al.\ 


2006|;|Evans|||2007|l 


|Kemp et a«.|||2010|l; 


Klimo et aZ. |(|2007|l; 



and references therein. 




FIG. 7 (Color online) Cartoon showing some of the possible 
acceleration mechanisms in the interaction with a thick solid 



target, including TNSA at the rear side (Sec II.C.l I, hole bor 



ing RPA at the front side (Sec II.C.2I, and backward accele r- 
ation in the plasma blow-off (see e.g. Clark et al. ( 2000b l). 
Also shown are the hot electron flow leading to sheath forma- 
tion and expansion at the rear side, and the associated return 
current. 



surface. Since Eg must backhold electrons with a typical 
"temperature" T/j, the typical spatial extension of the 
sheath Lg will be related to Eg by 



eEg 



Th 
Lg' 



(13) 



From dimensional arguments, assuming a steep inter- 
face and Tih and Th as the only parameters, Lg may 
be roughly estimated as the Debye length of hot elec- 
trons, Lg ^ Xoh — (?)i/47re^n/i)^/^. Assuming the sim- 
ple scalings of Sec |II.B| for Th, taking a laser irradiance 
IX^ = 10'^" W cm^ and a fractional absorption rjh = 0.1 
we find Th ~ S.lweC^ = 2.6 MeV, n,, - 8 x 10^° cm'^, 
Xdh = 4.2 X 10-5 cm and - 6 x 10^° V cm'^ This 
huge field will backhold most of the escaping electrons, 
ionize atoms at the rear surface and start to accelerate 
ions. As a rough estimate, a test ion crossing the sheath 
would acquire the energy Si ~ ZeEgLg — ZTh, resulting 
in MeV energies and a scaling as /^/^ if Th — £p given 
by Eq.([6]) holds. Protons from a thin layer of hydrogen- 
containing impurities on the surface will be in a very 
favorable condition for acceleration because of both their 
initial position, located at the maximum of the field, and 
their highest charge-to-mass ratio so that they will be 
more rapid than heavier ions in following electrons and 
screening the sheath field. This is the qualitative sce- 
nario for TNSA of protons as introduced by [Wilks et al.\ 
(2001) to explain their experimental results on proton 



acceleration (Hatchett et al. 2000 Snavely et al. 20001. 



The essential features of the TNSA mechanism have 
been supported by several experiments and TNSA has 
become the reference framework to interpret observations 
of multi-MeV protons from the target rear side. Various 
schemes for beam optimization and control have been 
designed on the basis of TNSA. A detailed discussion 



10 



of main experimental findings is reported in Sees III. A 
UlLBl and HTLEl 

From a theoretical viewpoint, there are two main cat- 
egories of models which describe TNSA, namely "static" 
and "dynamic" models which, depending on the start- 
ing assuptions, allow to provide simplified analytical de- 
scriptions to be used to interpret experimental data. 
These models and related numerical investigations are 



presented in Sec III.C 



2. Front surface acceleration 

Already in the first measurements of proton accelera- 
tion in the forward direction, the possibility of a contribu- 
tion originating at the front surface of the target was also 



conceived ( Clark et al. 2000a Maksimchuk et al. 2000 1 . 



As a consequence, mechanisms leading to ion acceleration 
in such region have also been extensively investigated. 

At the front surface, the intense radiation pressure of 
the laser pulse pushes an overdense target inwards, steep- 
ening the density profile and bending the surface; this 
process is commonly named as "hole boring" . The reces- 
sion velocity vub of the plasma surface may be estimated 
by balancing the electromagnetic and mass momentum 
flows, I/c ^ ni{miVhb)vhb- This corresponds to an en- 
ergy per nucleon Si = ^ I/{Anic). The scaling 
with the laser intensity / is more favorable than the /^/^ 
scaling for TNSA, and suggests Radiation Pressure Ac- 
celeration (RPA) effects to become more important for 
higher intensities. More accurate, relativistic and dy- 
namic modeling is presented in Sec |IV.A.l| along with 
related experimental indications. 

Radiation pressure action and hot electron temper- 
ature may also lead to the generation of coUisionless 



shock waves (Tidman and Krall 19711 with high Mach 
number M. Such waves are associated to the reflec- 
tion of ions from the shock front, resulting in a velocity 
Vi = 2Mcs and an energy per nucleon £i = 2mpM'^c1 = 
2{Z/A)M^Th, being = y^ZTijArrip the ion sound ve- 
locity. Such CoUisionless Shock Acceleration (CSA) sce- 
nario and related experiments are discussed in Sec |IV.B 



Finally, the possibility of front side (or bulk) acceler- 
ation being favored by resistivity effects is discussed in 



3. Acceleration in ultrathin, mass-limited and low-density 
targets 

Both TNSA and RPA may have different features in 
targets having peculiar geometrical and physical proper- 
ties, if compared to the sohd targets used in the 2000 ex- 



periments ( Clark et al. 2000a Maksimchuk et al. 2000 



Suavely et al. 2000 1 which were several micron thick and 



investigations of "ultrathin" , sub-micrometric targets re- 
quires extremely "clean", prepulse-free pulses to avoid 
early target evaporation and thus became possibile only 
recently thanks to the development of advanced tech- 



niques (see Sec III.E). The use of "mass-limited" targets 
which also have limited lateral dimensions (in the sub 
mm-range) allows the refluxing and concentration of hot 
electrons in a small volume, and may lead to higher ion 
energies via TNSA. These studies will be presented in 
SecHEEl 

For RPA, a sufficiently thin foil target is expected to 
be accelerated as a whole. Assuming the foil to be a 
perfect mirror of thickness £, its non-relativistic motion 
may be simply described by the equation miTiiidV / dt = 
21 /c from which we obtain an energy £i — mpV'^/2 — 
{2/mi){F/£c)'^ where F = J Idt is the laser pulse fluence. 
This is the basics of the "Light Sail" (LS) regime of RPA 



(Sec IV. A. 2 1 that seems very promising for the foreseen 



fast scaling and the intrinsic monoenergeticity. 

For extremely thin (a few nm) targets, the break- 
through of the laser pulse through the foil due to rela- 
tivistic transparency may stop LS-RPA, but at the same 
time lead to strong heating of electrons. This effect opens 
up a regime of enhanced acceleration possible, which has 
been also named Break-Out Afterburner (BOA) and will 
be discussed in Sec lIV.Cl 

In general, reducing the effective size of the target al- 
lows for laser pulse penetration, volumetric heating, and 
energy confinement, possibly allowing for efhcient ion ac- 
celeration even at low laser pulse energies. As a famous 
example, the interaction of ultrashort, moderate inten- 
sity (~ 10^^ W cm^^) pulses with sub- wavelength clus- 
ters allowed acceleration of ions up to energies sufficient 



to produce nuclear fusion reactions (Ditmire et al. 1997 



1999). A limitation on the use of such clusters as ion 



sources is the isotropic ion emission and the resulting 
low brilliance. "Droplet" targets with size of the order 
of one wavelength have been investigated as a trade-off 
approach, as discussed in Sec |III.E"!2l 

As mentioned above, special target materials may 
be used to produce plasmas with density close to ric 
(for laser wavelengths A ^ 1 ^va) in order to enhance 
the generation of hot electrons which drive TNSA (see 
SecIV.Dl. Gas jet targets have been also used both 
with A 1 /im lasers to investigate ion acceleration 



much wider than the laser spot diameter. Experimental 



in underdense plasmas (Sec |IV.D ) and with CO2 lasers 
(A ~ 10 /J,m) for studies of RPA and CSA in moderately 
overdense plasmas (see Sees IV. A and IV. B). Apart from 
the possibility to vary the background density, using flow- 
ing gas jets as targets is of interest because they enable 
the interaction with a pure proton plasma and are suit- 
able for high repetition rate operation as needed for most 
foreseen applications (Seclvl). 
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D. Particle-in-cell simulations 



The particle-in-ceU (PIC) method (Birdsall and Lang- 



don 1991 Dawson[ 1983), already mentioned in the In- 
troduction, is the most widely used approach to the 
kinetic simulation of plasmas. The PIC method pro- 
vides a solution to the Maxwell- Vlasov system using a 
Lagrangian approach, with fields and currents allocated 
on a fixed grid and the phase space represented by an 
ensemble of computational particles. Thus, the PIC 
method is mostly appropriate to describe collisionless 
laser-plasma interaction dynamics, although models are 
available either to implement either collisions [see e.g. 



Fi- 



uza et al. (2011 1 and references therein] or ionization [see 



e.g. Petrov et al. (20091 and references therein]. 

PIC simulations of laser interaction with solid density 
plasmas at peak densities typically exceeding lO^ric are 
a very demanding task even when most powerful super- 
computers are used. At least, one has to resolve temporal 



scales ■ 



and spatial scales ~ c/ujp where ujp 



-1/2 



is the plasma frequency. Thus, when approaching pa- 
rameters of a real experiment, relevant lenghts such as 
the laser beam waist may correspond to thousand of 
gridpoints in each spatial direction, and typical dynamic 
times to thousands of timesteps. In addition, kinetic ef- 
fects such as generation of hot electron tails in the dis- 
tribution function and large density variations need very 
large numbers of particles to be properly resolved. For 
these reasons, "realistic" 3D simulations with proper res- 
olution are typically beyond computational possibilities. 
This issue forces most of the times either to use a re- 
duced dimensionality or to relax the actual parameters 
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FIG. 8 (Color online) Results from multiparametric 2D simu- 



2006 



20021. 



lations for a double layer target ( Esirkepov et al. 
The maximum energy of protons accelerated from the rear 
layer is shown as a function of laser pulse energy £l and for 
different values of the intensity I, pulse length Lp and fo- 
cal spot diameter D. The target density and thickness are 
Ue = lOOric and i = X, respectively. Reprinted figure with 
permission from [Esirkepov et al. ( |2006| , Phys.Rev. Lett. 96, 
105001. Copyright 2013 by the American Physical Society. 



to some extent, e.g. by assuming relatively low densities 
or short scales. For some peculiar problems, development 
of hybrid modeling may be appropriate, as discussed in 
Sec lIII.C.4l 

Despite the above mentioned limitations, several 
groups have been able to perform large scale 3D simu- 
lations relevant to ion acceleration regimes such as, e.g.. 



TNSA (Pukhov 20011, RPA (Esirkepov et al. 2004 



Tamburini et aL[ |2012] and BOA ( |Yin et al.\ |201l"af . 



Use of parallel supercomputers has also allowed extended 
multi-parametric studies (Figjsj) in order to infer scaling 
laws and to evidentiate a transition from TNSA to RPA 



dominance at high intensities ( jEsirkepov et al. 20061. 
These and other numerical results will be discussed in 
SecsHniandllVl 



E. Ion diagnostics 

The specific properties of laser-driven ion beams (e.g. 
broad spectrum, high flux, significant divergence) have 
required either modifications of established diagnostics 
techniques or the development of new ones. 
Radiochromic film (RCF) ( [McLaughlin et "al. 



1996) 



is a detector which is favoured by many experimenters. 



since the early work by Suavely et al. (20001, mainly 



due to simplicity of use. This is a high-dose, high- 
dynamic range film, widely used in medical context for 



X-ray dosimetry ( |Niroomand-Rad et al. 1998). The films 



consist of one or more active layers containing a micro- 
crystalline monomeric dispersion buried in a clear plastic 
substrate. Typical examples are the IID810, MD55 and 
EBT2 Gafchromic varieties. After interaction with ioniz- 
ing radiation, the active material undergoes polymeriza- 
tion and the film changes its colour from nearly transpar- 
ent to blue. The consequent change in Optical Density 
can be calibrated against the dose released in the film, 
and therefore provide information on the fiux of particles 
directed at the layer. Tipically, RCFs are used in a stack 
arrangement, so that each layer acts as a filter for the 
following ones in the stack. The signal in a given layer 
will be due only to ions having energy E > Eb, where 
Eb is the energy of the ions which reach their Bragg peak 
within the layer (see Figj9]). In first approximation, for 
an exponential-like spectrum as those typically produced 
by TNSA, the dose deposited in a layer can be taken as 
proportional to the number of protons with E ~ Eb, 
allowing a rough spectral characterization of the beam. 
Various, more refined procedures have been developed for 
deconvolving the spectral information (either integrated 
across the beam or angularly resolved) in multilayer RCF 



(|Breschi et al. 


2004| Hey et al. |2008 


Kirby et al. 


Nuernberg et al. 




2009|). 



Plastic track detectors such as CR39 (Enge 19951, 
which have the advantage of being insensitive to X-rays 
and electrons, have been widely used also in multi-stack 
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FIG. 9 (Color online) A typical RCF stack obtained in an 
experiment with the TARANIS laser at Queen's University, 
Belfast. Higher numbers correspond to higher energy pro- 
tons with layer 1 corresponding to ~ 1 MeV and layer 9 to ~ 
12 MeV. Below each film the corresponding Bragg peak en- 
ergy for prot ons is reported. Reprinted figure from |Dzelzainis| 
et al. (20101, Las. Part. Beams 28, 451. By permission from 



Cambridge University Press (2013). 



arrangements, e.g. by Clark et al. (2000a). CR39 lay- 
ers (typically 0.25 — 1 mm in thickness) require etching 
in a NaOH solution after exposure to ions, so that the 
damage tracks created by the particles can be revealed 
thanks to the different etching rate in the track compared 



to the undamaged bulk ( Seguin et al. 2003 1 . After etch- 
ing, the single tracks can be counted, which provides a 
direct measurements of the number of protons hitting the 
detector. CR39 works better for low particle fluxes, as 
at high flux (tipically above ^ 10^ particles cm~^) or for 
long etching times the tracks start to overlap, leading to 



saturation ( ^Gaillard et al. 20071 



The interaction of laser-driven high-energy ions with 
secondary targets can initiate a number of nuclear re- 



actions (Sec V.El, which can been used to diagnose the 
beam properties with the ability to provide absolute par- 
ticle numbers with a linear response and virtually no sat- 
uration at high flux. The ^'^Cu(p, nf^Tni reaction in cop- 
per stacks has been used to quantify the proton numbers 
through measurement of j3~^ decay of ^'^Zn nuclei, using a 



Nal detector-based coincidence counting system ( Santala 



et al. 2001 Spencer ei aL 2001 ). Techniques employing 



a single Cu layer, in which a range of isotopes resulting 
from proton-induced nuclear reactions is analysed in or- 



der to reconstruct the proton spectrum, have also been 
used (Yang ei a/. 2004b I. Spectral resolution is provided 



by a combination of filtering and known thresholds for 
the considered reaction. The above described approaches 
provide particle flux integrated over the whole beam cross 
section. Contact radiography of (p, ?i)-generated isotopes 
in activation samples (where the activated foil is placed 



2008 



in contact with RCF) has been developed ( Clarke et al. 
Offermann et al\ [2010 ) as a way to achieve 2D 



images with high spatial resolution and extremely high 
dynamic range. Neutron spectra produced through fu- 
sion reactions of the type D(d, n)^He have been used as 
a diagnostic of laser-driven deuterium ions inside a laser- 
irradiated target (Habara et al. 2003 2004a|b ). 

Obtaining spectra with high energy resolution requires 
the use of magnetic dispersion techniques. In simple mag- 

the 



Chen et al. (2008a 



netic spectrometers [see e.g. 
ions, spatially selected by an entrance slit or pinhole, are 
dispersed along one spatial direction according to their 
energy by a ~ 1 T magnetic field B This arrangement, 
which discriminates particles according to their energy 
but not to their charge/mass ratio, is adequate for di- 
agnosing the high-energy proton spectrum in "standard" 
TNSA experiments in which protons are the dominant 



2002) 



accelerated species (Hegelich et al. 

A more complete spectral characterization of multi- 
component ion beams can be obtained with Thom- 
son Parabola spectrometers, based on the principle for 
mass spectrometry introduced by [Thomson ( |1911 ). A 
schematic of the device is shown in Fig|10[a). Ions are 
deflected by parallel E and B fields (with S ~ lO'' V/m 
(Sakabe et al. 1980)) resulting in a characteristic deflec- 
tion pattern in which species with different charge/mass 
ratio form separated parabolic traces in the detector 
plane, as shown in the typical image of FigjlOjb). Modi- 
fied magnetic and Thomson spectrometers, having imag- 
ing and angular resolution capability, have been also 



developed (Chen et al. 


2010 


Jung et al. 


2011a 


Ter- 


Avetisyan et al. 


20091. 



The detectors used in conjunction with these spectrom- 
eters are typically either CR39, Image Plates (IP), scin- 




FIG. 10 (Color online) Left frame: schematic of a Thom- 
son parabola (courtesy of S. ter- Avetisyan) . Right frame: 
a typical example of ion traces obtained with the Thomson 
parabola. 
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tillating plates or MicroChannel Plates (MCPs). Photo- 
stimulable IPs are film-like radiation image sensors, de- 
veloped for X-ray medical imaging, which are composed 
of specially designed phosphors that trap and store radia- 
tion energy in metastable excited states, and can be abso- 



ated, on a longer time scale, if the proton number is not 
high enough to balance the charge of the escaping hot 
electrons, and especially if impurity protons are removed 
before the interaction, for example by pre-heating the 
target (HegeHch et al. 2002 1. In this way, ions of sev- 



2011 Mancic et al. 2008). Scintillating plates (Green 2005) 



et al. 2010) or MCPs ( Ter-Avetisyan et al. 2005) are 



lutely calibrated in terms of particle flux ( Freeman et al. eral different species may be accelerated ( Hegelich et al. 



favoured in situations where online detection is required 
(e.g. high-repetition laser systems), as the scintillator 
screen or the MCPs phosphor are imaged on a CCD and 
the detector does not require replacing after exposure. 
Scintillators can also be used for beam profiling (Sa kakij 



et al. 20101, with potential for energy range selection 



(Green et al. 20111. 



A different approach also allowing online beam mon- 
itoring is the use of Time of Flight (TOF) techniques, 
where the broadband ions are left to propagate over a 
given distance and then detected employing scintillating 



plates coupled to a Photo Multiplier (Nakamura et al. 



2006), Faraday Cups or semiconductor detectors (Mar- 



garone et al. 2011). The time- varying signal produced 



by the detectors maps the ion energy spectrum, although 
the finite response time of the detector and realistic prop- 
agation distances limit the use of these techniques for 
measurements up to a few MeV/nucleon energy. State- 
of-the art TOF-MCP detectors allow for measurements 
of protons with a kinetic energy up to ~ 20 MeV/nucleon 



(Fukuda et al. 2009) 



III. TARGET NORMAL SHEATH ACCELERATION 



A. TNSA Scenario. Main experimental observations 



As anticipated in Sec II. C the TNSA process (Wilks 



et al. 2001 1 is a consequence of the huge charge sepa- 



ration generated by hot electrons reaching the rear side 
of the target. There, a cloud of relativistic electrons is 
formed, extending out of the target for several Debye 
lengths, and giving rise to an extremely intense electric 
field, mostly directed along the normal to the surface. A 
consequent distinctive feature is that ions are accelerated 
perpendicularly to the surface, with high beam collima- 
tion. The electric field generated at the rear surface will 
depend on parameters of the electron distribution (tem- 
perature, number, divergence) as well as parameters of 
the surface itself (mostly its density profile, as detailed 
below). 

The acceleration is most effective on protons, which 
can be present either in the form of surface contaminants 
or among the constituents of the solid target, as in plastic 
targets. The heaviest ion populations provide a positive 
charge with much more inertia, thus creating the charge 
separation which generates the accelerating field. Part 
of the heavy population can be also effectively acceler- 



Several observations strongly supported the TNSA sce- 
nario taking place at the rear side. Already |Snavely et al. 
( 2000 ) gave clear evidence that the emission was normal 



to the rear surface using wedge targets which effectively 
have more than one rear surface. Two separate proton 
beams were observed in the directions normal to the two 



rear surfaces of the wedge (Fig 111. 

Mackinnon et al. ( 2001[ ) reported experimental obser- 



vations of the interaction of ultraintense laser pulses us- 
ing targets with and without preformed plasmas on the 
rear surface of the foil. The peak and mean energies of 
the proton beam were found to depend strongly on the 
plasma scalelength at the rear of the target. While an en- 
ergetic proton beam was obtained with an unperturbed 
rear surface, no evidence of high energy protons was 
recorded when a large local scalelength in the ion density 
at the rear surface was induced, consistently with the de- 
pendence of the accelerating field on the scalelength in 



Eq.(13) 



Hegelich et al. (20021 used as targets Al and W foils. 



resistively heated to remove Hydrogen contaminants and 
coated on the rear side with thin C and CaF2 layers, re- 
spectively. The observation of high-energy C, Ca and F 
ions out of these prepared source layers prove the exis- 
tence of an effective rear-surface acceleration mechanism 



(Fig 12 1 . Further evidence was given by the work of Allen 



et al. (2004) who showed that, removing contamination 



from the back surface of Au foils strongly reduced the 
total yield of accelerated protons, while removing con- 




10G 



10 



Q.1 
06 



FIG. 11 (Color online) Proton emission from a wedge target 
effectively having two rear surfaces. Two separate spots are 
produced on the detector, showing that most of the protons 
originate from the rear side of the target. Reprinted figure 
with permission from Snavely et al. (20001, Phys. Rev. Lett. 
85, 2945. Copyright 2013 by the American Physical Society. 
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tamination from the front surface of the target had no 
observable effect on the proton beam. 

A further proof that ions are accelerated at the rear 
side was given by the observation that a structuring (i.e. 
grooving) of the rear surface produced modulations in 



the proton beam (Cowan et al. 2004). This effect also 



evidences the high laminarity of the beam and allows to 
measure its emittance, as discussed in Sec jllLBl 

Direct experimental evidence of the generation of an 
initial intense sheath field at the rear surface and a late 
time field peaking at the beam front was provided by'Ro-l 



magnani et al. ( 2005 ) using the proton imaging technique 



(see Sec V.Al. In other words, TNSA provided itself an 
unique diagnostic which allowed a direct experimental 
confirmation of the essential nature of the acceleration 



process. Fig 13 shows a temporal series of "proton im- 
ages" in which the propagation of the bell-shaped front 
of ion expansion can be observed. 

More recent developments have been obtained playing 
with the detailed properties of both the laser pulse and 
the solid target, showing the possibility to significantly 
control and optimize the TNSA process, and the capa- 
bility of achieving interesting and promising variations 
of the main scheme, also in the light of possible specific 



applications. These topics will be presented in Sec III.E 
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FIG. 12 (Color online) Effect of impurity removal on Carbon 
ion spectra. Frames a) and b) show C ions traces (from CR- 
39 tracke detectors) and spectra from Al foils coated with a C 
layer on the rear side, in the presence of hydrocarbon contam- 
inants on the surface. In frames c) and d), the contaminants 
had been previously removed by resistive heating. Reprinted 
figure with permission from |Hegelich et aZ.| ( |2002| ), Phys. Rev. 
Lett. 89, 085002. Copyright 2013 by the American Physical 
Society. 



B. Beam properties characterization 

Several experiments have investigated in detail the 
properties of the TNSA ion beams. The energy spec- 
trum of the beams is typically broadband, up to a cut-off 
energy (see Fig|2|. The particle number per MeV can 
be roughly approximated by a quasi-thermal distribution 



with a sharp cutoff at a maximum energy (Fuchs et al. 



2005 Kaluza et al. 2004) which scales with the laser 



parameters as will be discussed in Sec. |III.D| Many ex- 
periments have reported spectral observations for a wide 
range of laser and target parameters, and partial surveys 



have been provided in a number of publications ( Borghesi 



et al. 2006 Fuchs ei aL 2006b Zeil et al. 2010). 



A number of experimental studies have been devoted 
to the investigation of the spatial and angular character- 
istics of the emitted beams. These latter are closely de- 
pendent on the electron sheath spatial distribution, and 
consequently on the target properties (resistivity, surface 
roughness, etc) affecting the electron propagation. 

It was observed relatively early on that the use of con- 
ducting targets leads to smooth proton beam profiles 
with a sharp boundary, as detectable, for example, in 



RCF images (Fuchs et al. 2003 Snavely et al. 20001, 



while using dielectric targets creates non-homogeneities 
in the proton density across the beam section. In the lat- 
ter case, the transport of the electron current is prone to 
electromagnetic instabilities, which break the hot elec- 
tron flow into filaments. This leads to an uneven elec- 



tron sheath at the target rear (Manclossi et al. 2006) 



and consequently to a modulated proton beam cross sec- 
tion (iFuchs et al.\ 120031 iRoth et al.\ 120021). The close 



correlation between proton beam properties and electron 
beam transport characteristics has indeed been exploited 
in a number of experiments, which have used the proton 
beam as a diagnostic for the electron beam behaviour in- 
side the target, revealing, beside the aforementioned dif- 



ferences related to the target conductivity (Fuchs et al. 



2003 ) , effects of magnetic coUimation on the beam trans- 



port (Gizzi et al. 2011 Yuan et al. 2010) or the role 



of lattice structure in dielectric targets ( McKenna et al. 



2011). Other factors that can lead to structured beam 



profiles even in conducting targets are surface roughness 
at the target rear, resulting in a randomized local orien- 
tation of the protons ( Roth et al. 2002 1 , and intensity 



modulations in the focal spot which can be coupled to 
the protons via structured electron beams in medium-Z 



thin targets ( Fuchs et al. 2003 1 



The existence of a sharp angular boundary in the pro- 
ton angular distribution (clearer in higher-Z and thicker 
targets) is consistent with a bell-shaped transverse distri- 
bution of hot electrons in the rear surface sheath due to 
the fact that the density will naturally be higher along the 
laser axis and decrease with transverse radius. Protons 
are accelerated normal to the local iso-density contour, 
and the presence of an inflexion point in the sheath there- 
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FIG. 13 Proton probing of the expanding sheath at the rear surface of a laser-irradiated target. Frame a): set-up fo the 
experiment. A proton beam is used as a transverse probe of the sheath. Frames b)-g): Temporal series of images produced 
by to the deflection of probe protons in the fields, in a time-of-flight arrangement. The probing times are relative to the peak 
of the interaction. Frame h): a deflectometry image where a mesh is posed between the probe and the sheath plasma for a 
quantitative measure of proton deflections. Reprinted figure with permission from Romagnani et al. 
95, 195001. Copyright 2013 by the American Physical Society. 



(20051, Phys. Rev. Lett. 



fore results in a maximum angle of acceleration (Fuchs duced by Cowan et al. (2004), has successfully been used 



et al. 2003). Comparison of experimental data with sim- in several experiments for diagnosing the emission prop- 



pie electrostatic models indicate that the shape of the erties of the beam ( Nuernberg et al. 2009 1 . From these 



accelerating sheath is generally Gaussian (Carroll et al. 



2007 Fuchs et al. 2003 ) as also observed directly in 



sheath imaging data (Romagnani et al. 2005 ), see Fig 13 



A modulation of the proton beam angular distribution 
can be introduced purposefully by nanostructuring the 
target surface. A technique based on micro-machining 
shallow grooves on the rear surface of the target, intro- 
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FIG. 14 Top frames: modulations in the proton distribu- 
tion, for different energies, on a RCF detector from a target 
with micro-grooves imprinted on the rear side. The target is 
a 18 fim thick Al foil irradiated at 10^^ W cm~^. Bottom 
frames: effect of electron removal by magnetic fields, showing 
that the proton beam emittance is not significantly affected. 
The images are for 6.5 MeV protons and the target thickness 
is 40 ^m. Reprinted figure with permission from |Cowan et aL\ 
( 2004[ ), Phys. Rev. Lett. 92, 204801. Copyright 2013 by the 
American Physical Society. 



patterned targets, a periodic modulation of the beam an- 
gular envelope arises during TNSA due to the local per- 
turbation of the target normal direction, which causes 
an initial beam microfocusing at the groove locations. 
As the sheath expands, the local modulations are added 



over the global divergence of the beam ( Ruhl et al. 2004 ) 



and are observable as a modulation of the proton dose on 



the detector (Fig 14). The modulations can be used as a 



spatial fiducial from which one can infer the dimensions 
of the area from where ions are accelerated, i.e. the pro- 
ton or ion source size ( Brambrink et al. 2006 ) . Similar 



information has been obtained by considerations based 
on the projection by the ion beam of patterned objects 



(Schreiber et al. 2004). 



e.g. metal meshes (Borghesi et al.. 2004) or knife-edges 



A crucial property of laser-driven ion beams is their 
laminarity. In an ideal laminar source, there is a linear 
correlation between the radius within the source from 
where a particle is emitted and its angle of emission. 
The degree of laminarity of charged-particle beams is 
typically expressed in terms of their transverse emit- 
tance, a quantity which is proportional to the area of 
the bounding ellipsoid of the distribution of particles in 



phase space (Humphries 1990 1. The highest quality ion 



beams have the lowest values of transverse and longitu- 
dinal emittance, indicating a low effective transverse ion 
temperature and a high degree of angle-space and time- 
energy correlation, respectively. Transverse emittance 
has been measured in a number of experiments. Meth- 
ods based on mesh projection (which is broadly equiv- 
alent to the established "pepper-pot" method) indicate 



that e < O.Itt mm mrad (Borghesi et al. 2004 Ceccotti 



et al. 2008 Nishiuchi et al. 20081. The above discussed 



groove imaging technique allows a full reconstruction of 
the transverse phase space, and possibly a more pre- 
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cise estimation of the transverse emittance ( 


Brambrink 


et al. 


2006 


Cowan et al. 


2004 


Nuernberg et al. 


2009) 



which, for protons of up to 10 MeV, has been estimated 
as 0.004 mm mrad, i.e. 100-fold better than typical RF 
accelerators and at a substantially higher ion current (kA 
range) . 

It has also been found that the removal of the co- 
moving electrons after 1 cm of the quasi-neutral beam 
expansion did not significantly increase the measured 



proton transverse emittance, as shown in Fig 14 (Cowan 



et al. 2004 1 . This last observation is important since, in 
order to take advantage of the exceptionally small proton 
beam emittance in future applications, e.g. to capture 
them into a post-accelerator, removal of the co-moving 
electrons without significantly perturbing the protons is 
crucial. 

The ultra-low emittance stems from the extremely 
strong, transient acceleration that takes place from a 
cold, initially unperturbed surface and from the fact that 
during much of the acceleration the proton space charge 
is neutralized by the co-moving hot electrons. Using the 
ion beam as a projection source, having a low-emittance 
beam is equivalent to projecting from a virtual point-like 
source located in front of the target, with much smaller 
transverse extent than the ion-emitting region on the tar- 
get surface (Borghesi et al. 2004 Nuernberg et al. 2009). 
As will be discussed in Sec V.A[ this property of laser- 



driven ion beams allows to implement point-projection 
radiography with high spatial resolution. 



C. TNSA modeling 

The experimental observations and the considerations 
summarized in Sec |III.A| suggest the following starting 
assumptions, leading to the formulation of a relatively 
simple system of equations which can be investigated an- 
alytically and numerically (Passoni et al. 2004). First 



of all, we assume an electrostatic approximation, so that 
the electric field E — —'Vcj) where the potential (f> satisfies 
Poisson's equation 

V^(/) = 47re(ne — ZjUj), 



(14) 



with the sum running over each species of ions, hav- 
ing density Uj and charge Zj . As a consequence of the 
laser-solid interaction, the electron density rig may be 
described as composed of at least two qualitatively dis- 
tinct populations, which will be labeled cold and hot 
in the following, having densities Uc and rih such that 
He = ric + Uh. In the simplest approach, thermal effects 
are neglected for the cold population, while is given 
by a one-temperature Boltzmann distribution, 

rih = no/.e^*/^'^ . (15) 

This expression can be a reasonable first approximation 
to account for the presence of the self-consistent sheath 



field and has actually been used in many works on TNSA[^ 
but, as discussed below, it can lead to serious problems 
when the main goal is the estimation of the maximum 
energy of the accelerated ions. Alternatively, the elec- 
tron dynamics can be included via either fluid or kinetic 
equations. It is most of times appropriate to consider 
two different ion species, a light (L) and a heavy {H) 
population; in this way it is possible, e.g. to model the 
acceleration of light species present on the surfaces of a 
solid target made of heavy ions. 

Depending on the description of the ion populations, 
two main categories of TNSA models, to be discussed in 
detail in Sees. III.C.l and III.C.2 respectively, may be 



identified as follows. The first includes static models in 
which it is assumed that the light ions, or at least the 
most energetic ones, are accelerated in the early stage of 
the formation of the sheath, so that the latter may be 
assumed as stationary. In these conditions, the effects of 
the light ions on the electrostatic potential are usually 
neglected, while the heavy ion population of the target 
is considered immobile. The aim is thus to provide the 
most accurate description of the sheath depending on as- 
sumptions on the hot electron distributions. The second 
category includes dynamic models where the system is 
described as a neutral plasma in which the ions acquire 
kinetic energy in the course of the sheath evolution. In 
several cases a unique ion component is considered. This 
approach is strongly connected to the classic problem of 



plasma expansion in vacuum, first considered by Cure- 



vich et al. (1966). In a cold fiuid description, neglecting 
relativistic effects the ions are described by the equations 

dt 
dnj_ 

at 



V-(7ijU,) = {]=L,H), 



(16) 
(17) 



where = Uj(r,t) is the fluid velocity. If the ions are 
described kinetically, their Vlasov equation for the phase 
space distribution fj = fj{r,-v,t) is 



dt 



V/, 



Z,e„ 



(9v 



0. 



(18) 



Most of general studies of plasma expansion and related 
ion acceleration developed both before and after TNSA 
experiment as well as more specific models of TNSAj^ 
so far proposed in the literature can be considered as 



See e.g . [Albright et ~ai\ \200Q\ ; |Mora| (|2003|l; 'Nishiuch i et al. 



Robinson et al. (2006). 



Passoni et aZ.| ( |2004| ; and' 

st udies o t plasma expansion m vacuum mcludc Allen and 

Andrews! (|1970|l; |Cro w a/.| (|T975|; |Denagt | (1979 ); Do rozh~ 
ina and Semcnov (1998); Gurevich et al. (1966); Mora and 
PcUat (19791; Pcarlman a nd Morse, ( ,1978,) ; and W 
(1971). Early papers focused on modeling ion acceleration in 
laser-produced plasmas include [Kishimoto et al\ ( |T983| ) ; |Peari^ 
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suitable simplifications of the previous equations, falling 



physically motivated assumptions. Most of these models 



approximation, and planar in most cases. This latter as- 



sumption, when applied to TNSA modeling, requires the 



rear surface to be sufficiently flat and the electron cloud 



to be spatially uniform in the plane normal to the ion 



into one of the two above mentioned categories (or suit- 



able combinations of them) and obtained adding further. 



assume a ID geometry, consistently with the electrostatic 



motion- 



Notice that all the models proposed to describe TNSA 



are, to a large extent, phenomenological, i.e. they need as 



input parameters physical quantities which are not pre- 



cisely known. Since these descriptions give a simplified 



picture of the acceleration process, the "best" model in 



this context may be considered the one which provides 
the best fit of experimental data with the lowest set of 



laser and target parameters. This issue will be discussed 



in SecjllLD In principle, these difficulties could be over- 



come performing "realistic" numerical simulations, but 
actually also the latter always consider "model" problems 
because of intrinsic difficulties in the numerical study of 



these phenomena, such as for example the large varia- 



tions of density from the solid target to the strongly rar- 
efied expansion front. At present, a complementary use 



of simple models, presented in Secs|III.C.l||III.C.3[ and 



advanced simulations, discussed in Sec III.C.4 seems the 



most suitable option to theoretically approach TNSA 



|1. Quasi-static models 



suming, on the time scale of interest (i.e. in the sub- 



ps regime), immobile heavy ions, an isothermal laser- 



produced hot electron population, and the light ions to 



Starting approximations of static models consist in as- 



be sufficiently few to neglect their effect on the evolu- 



tion of the potential so that they can be treated as test 



particles. In this limit, if Eq.(15l is used to describe hot 



electrons and neglecting thermal effects for cold electrons. 



the potential in planar geometry is determined by 





= 47renoh 




(19) 



where we assumed the background charge to fill the x < 



region with uniform density. The corresponding electron 



I man and Morse| | |1978[ l; |True et al.\ | |1981| l; and |Wickens et al.\ 
| |1978 I . More recent works stimulated by the TNSA experiments 
include Betti et aL|l| 2005');'Cecche rini et aZ.| ||2006|l; [Kovalev and| 

[ Bychenkov| ( |2003 ^;|Mora ( 2003^ , 2005^ ; and |Peano ei aL| | |2007[ |. 
A list ot papers dcscribing 'i'JNyA moc els mostly based on a static 
modeli ng include Albright et o Tj | j2006jl; [Lontano and Passoni| 
l|2006|l;|Passoni and Lont ano (2004 2008) ' Passon i et gf.] | |2004l l7 
[Robinson et aL| l |2006[ l; and |Schreiber et af.| l |2006( |. 



density and electric field can be calculated, as well as 
the energies of test ions moving in such potential. This 
can be considered the simplest self-consistent approach to 
theoretically describe the TNSA accelerating field. The 
solution of Eq.(|19|, in the semi-infinite region a; > is 



(Crow et al. 19751 



e 



In 1 



2eA 



Dh 



(20) 



where Xoh — \/Th/{4:ne^noh)- The field reaches its max- 
imum at the surface and is given by 



E{0) 



Eq . 



Eo 



Th 



eX 



Dh 



(21) 



which justifies the simple estimates used in Sec |II.C| 
However, the electrostatic potential (20) leads to an in- 



finite acceleration of a test proton which is initially at 
zero energy in x = 0. The reason is that the apparently 
reasonable choice of the Boltzmann relation poses sever 



difficulties to the analysis (Passoni and Lontano 2004) 



because, in order to have an electron density equal to 
zero at infinity, the self-consistent electrostatic potential 
must diverge at large distance from the target [mathe- 
matically, — ^ — oo as a; +oo, see Eq.([T5|]. This is 
not a pathological consequence of the one dimensional 
approximation but it is related instead to the fact that 
the Boltzmann relation implies the existence of particles 
with infinite kinetic energy, which is not physically mean- 
ingful [see also § 38 of Landau and Lifshitz (1980)]. This 
unphysical behavior can be avoided by assuming an up- 
per energy cut-off £c in the electron distribution function, 
so that e(/) — —£c as a; — > +oo and the electric field turns 
to zero at a finite distance. The cut-off assumption can 
be justified as a consequence of the laser-solid interac- 
tion producing electrons with a maximum kinetic energy 
and of the escape from the system of the most energetic 



ones (Lontano and Passoni 2006 Passoni and Lontano 



2008). Experimental indications of target charging due 



to electron escape have been found by, e.g. Kar et al. 



(2008b) and Quinn et al. (2009a). The finite range of 



the electric field driving TNSA is also apparent in direct 



measurements (Romagnani et al. 2005). 



As a first development, still using the Boltzmann rela- 



tion, it can be assumed the ID solution given by Eq.( 19 ) 
to hold only up to a longitudinal distance roughly equal 
to the transverse size of the sheath, because at larger dis- 
tances 3D effects should be taken into account, contribut- 

Al- 



2006) 



ing to remove the divergence ( Nishiuchi et al. 
ternatively, by assuming that the hot electron population 
occupies only a finite region of width h, the solution of 



Eq.(18) in the vacuum region < x < h, together with 



the corresponding electric field and electron density can 



be determined (Passoni and Lontano 2004). 



Another possibility, explored by Schreiber et al. ( 2006 ) 



has been to heuristically assume that the hot electrons 
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and b = Tc/Th, as shown in Fig|l5| ([Passoni et oH |2004 1 . 
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FIG. 15 Electric field profile in a sheath with two elec- 
tron temperatures. The field is normalized to Th/eXoh 
and is shown for cold-to-hot electron temperature ratio b = 
Tc/Th = 0.01 and for different values of the pressure ratio 
ab — poc/poh = 1 (dotted line), ab — 10 (dashed line) and 
ab = 100 (solid line). The x coordinate is normalized to the 
cold electron Debye length Xdc corresponding to ab = 10. 
Reprinted figure with permission from 



Phys. Rev. E 69, 
Physical Society. 



Passoni et al. ( 2004 \ 
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expansion in vacuum creates a cylindrical quasi-static 
cloud in the vacuum, behind the target, and a circular 
positive surface charge on its rear face. The generated 
electrostatic potential is evaluated on the symmetry axis, 
along which the most energetic ions are accelerated. The 
total surface charge and the radius of the distribution are 
model parameters estimated from experiments (see also 
SecinrDl). 



In order to consistently overcome the previously dis- 



cussed limits, Lontano and Passoni (2006) proposed to 



solve the Poisson equation by assuming that a quasi- 
stationary state is established where only those elec- 
trons {trapped electrons) with negative total energy W — 
mc^{'^ — 1) — 60 are retained, while those with positive 
total energy are lost from the system. The corresponding 
trapped electron density, given by rih — /ly^g fe{x,p)dp, 
is included in the Poisson equation and the correspond- 
ing analytical solutions determin ed ([Lontano and Pas 



soni 



|2006| [Passoni et al.\ |2010b[ [Passoni and Lontano| 



2008). As a general feature, the potential, the electro- 



static field and the electron density distributions go to 
zero at a finite position Xf oi the order of several hot 
Debye lengths. 

If both electron populations, hot and cold, are con- 
sidered, it is possible to find an implicit analytical so- 



lution of Eq.(19) both inside the target and in the vac- 
uum region. Using a two temperature Boltzmann re- 
lation to describe the electron density, that is — 
riQh exp {e(j)/Th)+nQc exp {e<j)/Tc), the electric field profile 
turns out to be governed by the parameters a = riQc/noh 



The presence of the cold electron population strongly 
affects the spatial profiles of the field, which drops al- 
most exponentially inside the target over a few cold elec- 
tron Debye lengths. An estimate of Tc, as determined 
by the ohmic heating produced by the return current 



(Sec II.B.3 ), is required. A simple analytical model of the 



process has been proposed (Davies 2003 Passoni et al. 



2004), to which we refer for further details and results. 



The quasi-static approach allows to draw several gen- 
eral properties of the accelerating TNSA field. The spa- 
tial profile is characterized by very steep gradients, with 
the field peaking at the target surface and decaying typi- 
cally over a few /im distance. The most energetic ions, ac- 
celerated in the region of maximum field, cross the sheath 
in a time shorter than the typical timescale for plasma 
expansion, electron cooling and sheath evolution. As a 
consequence the static approximation will be more ac- 
curate for the faster ions. Assuming a time-independent 
field also requires the electron cloud not to be affected 
by the ions flowing through it, which implies the number 
of the accelerated ions to be much smaller than that of 
the hot electrons, iVi <C -/Ve- A quasi-static model not re- 
quiring this assumption was proposed by [Albright et aL\ 
( 2006 ) who included effects of the accelerated ion charge 



on the electric field by modeling the layer of light ions 
(having areal charge density Ql) as a surface layer of 
density = {Q l / Z Le)S{x — x l) ■ Eq.(19l is then solved 



as a function of the instantaneous position x l ■ An exten- 
sion of this model, using an adiabatic descritpion of the 



hot electrons, was proposed by Andreev et al. (2008), to 



investigate the variation of the maximum light ion energy 
as a function of the heavy ion target thickness. 

On the basis of the above discussions we expect static 
models to be most reliable for the estimate the cut-off in 
the ion energy spectrum. This estimate requires as an 
input a few parameters, depending on the model. This 
issue will be discussed in Sec lIII.Dl 



2. Plasma expansion into vacuum 

A description of ion acceleration over long times 
and/or in conditions such that the quasi-static model- 
ing of Sec [III.CTT1 is not valid anymore demands for the 
inclusion of the ion dynamics. The description may be 



based either on a fluid modeling, using Eqs.(16 17), or 



on a kinetic one using Eq.(18) 



The simplest approach is obtained using a ID fluid 



approach, invoking quasi-neutrality, using Eq.(15) and 



assuming a single ion and electron population expand- 
ing in the semi- infinite space a; > 0. Eq.(19) is sub- 



stituted by the simpler condition rif, — ZiTii, the in- 
dex i denoting the single ion component. The boundary 
conditions are that the electron density should remain 
equal to the background value well inside the plasma, so 
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FIG. 16 Sketch of the density and velocity profiles from 
the self-similar solution for isothermal plasma expansion, 
Eqs.(22 13|. The front of charge separation at x = 
and the rarefaction front at a; = —Cst are also indicated. The 
electric field is uniform in the —Cgt < x < Xf{t) region. 



n^{—oo) = ng, and should vanish in vacuum far from the 
surface, ne{+oo) = 0. Together with Eqs. ( 16|17 1, the 
resulting system admits the classical self-similar solution 



first found by Gurevich et al. ( 1966 ) , 



no exp 



X 

cd 



1 



X 

t 



(22) 



where x/t is the self-similar variable, L = ni/\dxni\ = 
Cst is the local density scalelength, and the expressions 
are valid for x > ~Cst. The profiles corresponding to 



Eqs. (22 1 are sketched in Fig 16 



As a consequence of the quasi-neutral approximation, 
the physical quantities describing the plasma dynam- 
ics present several diverging behaviors, like the unlim- 
ited increase of Ui with x. This implies that the neu- 
tral solution must become invalid at some point, which 
can be estimated by equating the local density scale- 
length L to the local Debye length Ad- This provides 
Xf{t) = Cst[2\ii{ujpit) — 1], the corresponding velocity 
Uf — dxf/dt = 2csln{ujpit) and the electric field at 
the ion front, Ef = E{xf) — 2EQ/ujpit, where Eq = 
{4:TTnoThy^^ ■ This estimate gives twice the self-similar 
field E = Th/{eCst). The argument also defines the front 
of the fastest ions moving at velocity u / and thus it gives 
also the high-energy cut-off in the energy spectrum of the 
ions in this description. 

Eqs. (22) are also singular for i — ^ 0, i.e. at the earli- 
est instants of the expansion in which quasi-neutrality 
also breaks down. In general, in the sub-ps regime 
the inertia of ions is important and the assumption 
of quasi-neutrality must be consistently abandoned and 
ultimately a self-consistent analysis can be developed 



through numerical simulations (Sec III.C.4). Still assum- 
ing, for the sake of simplicity, that only a single ion popu- 
lation and a single-temperature Boltzmann electron pop- 
ulation are present, and ni{t = 0) = no8(— x), Eq.(20) 



can be used to define the initial conditions for the elec- 
tric field at the time t = aX which the ion acceleration 
process begins. The following interpolation formulas for 



the electric field and ion velocity at the ion front 



(23) 
(24) 

where r = ujpit/\^2e, give the correct behavior at t ^ 



e -I- 1 
Uf{t) ~ 2cs In (t + v/r^Tl) 



for both the electric field [see Eq.(21|] and the front 
velocity, and reduce to previous expressions for ujpit ^ 1. 



These formulas fit well numerical calculations by Mora 



(2003) using a Lagrangian fluid code. Related results of 
these studies using the fluid and the kinetic descriptions 
can be found in the literatureH 



The major drawback of Eq.(24) is that the maximum 
velocity of ions, and hence the cut-off energy, diverges 
logarithmically with time. This is not surprising, being 
an unavoidable consequence of the isothermal assump- 
tion and the chosen boundary conditions: the system has 
an infinite energy reservoir in the electron fluid and thus 
it is able to accelerate ions indefinitely. Nevertheless, the 



simplicity of Eq. (24 1 has proven to be attractive, thus it 
has been suggested to insert a phenomenological "max- 
imum acceleration time" tacc at which the acceleration 
should stop. Such a formula has been used in attempts 
to fit experimental data ( Fuchs et al. 2006b ) . We will 



III.D There is no eas- 



conie back to this point in Sec. 
ier way to remove this unphysical behavior from the ID 
planar model but to give a constraint of finite energy 
(per unit surface). In this way, the electron tempera- 
ture decays in time due to the plasma expansion and to 
coUisional and radiative losses. The electron cooling co- 
operates with the effects of finite acceleration length and 
maximum electron energy in the determination of a finite 
value for the maximum energy gain. 

The expansion of plasma slabs (foil) of finite thick- 
ness, and hence of finite energy, has been considered an- 
alytically and numerically. In these models The electron 
temperature is taken as a function of time, Th = Th(t), 
determined either by the energy conservation equations 



(Betti et al. 


2005 


Mora 


2005 




of Th{t) (Bychenkov et al. 


2004 



Analytical solutions for the plasma expansion can 
be found, in the quasi-neutral approximation, also for 



the kinetic Vlasov equation (18), using either the self- 



similar theory (Dorozhkina and Semenov 1998) or a re- 



normalized group theoretical approach (Kovalev et al. 



2001 2002). Two-temperature electron distributions 



have also been considered IHl 



See e.g. [Crow et aL|l|l975| ); [Denavit|(| l979');l Mora| | (2003l l; |Pearl-| 
I man and Mo rse| ||1978[l; and |Widner et af.,(,197l| f or the case of 
a single electron population a nd [Bychenkov et ar] (|2004); Mora] 
( |2005^ ; and ,Tikhonchuk et aZ.| | |2005[ l for "the case of two electron 

compon ents. 

13 See e.g. [Bezzerides et a;.| | |l978l l; |Diaw and Mora| pOllt ; |Gure-| 
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|3. Multispecies expansion 



] 



in which the dynamics of a heavy ion component is con- 



sidered in addition to light ions. The most pecuhar effect 
of the presence of two (or more) ion species, for appro 



priate parameters, is the appearance of spectral peaks, 



which are of interest both as a strong experimental sig- 



We now describe the expansion of a two species plasma. 



nature and for application purposes. 



The problem of two species expansion was studied by 



several authors since longtime ( Bezzerides et al. [1978 



mostly follow the more recent work of Tikhonchuk et al. 



Gurevich et al.\ 1973 Srivastava et al. 1988). Here we 



(2005), where a simphfied description is given based on 



the ordering assumptions 



Ah/Zh 



»1, 



TV : 



ZnnH 
Zlul 



> a. 



(25) 



small, and that the L ion plasma frequency is higher so 



that the dynamics of L ions is faster. These assumptions 



the H ions can be described as a single species expansion 



sound speed. L ions are treated as test particles in this 
region, where they are accelerated by the electric field 



E QC [see Eq.( 13[)], which is thus stronger than that 
would be created in the expansion of the L ions alone. 



The L ion velocity and density in this region can be ob- 
tained using the ID fluid, self-similar equations with the 



These conditions state that the H species is quite heav- 



ier than the L one, that the concentration of the latter is 



allow to assume that, near the rarefaction front, the ef- 



fect of L ions is unimportant and that the dynamics of 



as in Sec|III.C.2[ where the relevant parameter is the H 
ion sound speed cr = \J ZuThj Ahtti^ < cl, the L ions 



above given electric field ([Tikhonchuk et al. 2005|), ob 



taining for the velocity profile vl — Cl^/^ (1 -I- x/cHt) 



T7T 



Noticeably, the L ions velocity and density profiles vary 
slowly in space compared to the H ion ones, and the L 



ion flux is almost constant. Beyond the H ion front, only 
L ions are present and they can be described again by a 



single species expansion, vl — Cl + x/t [see Eq.(22)]. 



However, matching of the velocity profiles in the region 



behind the H ion front implies the existence of a transi- 
tion region where the velocity is approximately constant. 
This corresponds to a plateau region in the phase space 



and in a peak in the L ion energy spectrum. The heuristic 
reason for plateau formation is that the L ions are accel- 



erated more efficiently behind the H ion front than ahead 



of it. Fig 17] shows the velocity spectrum from numeri- 
cal results [Tikhonchuk et al. ( 2005 ) using a Boltzmann- 



Vlasov-Poisson model ( Bychenkov et al. 2004 ) based on 



Eqs. ( 14 ■ 15 ■ 18 ) , compared with analytical estimates from 



the self-similar solution. 




V/C. 



2 3 4 5 6 



FIG. 17 (Color online) Velocity spectrum of heavy ions (H, 
blue thick lines) and light ions (L, red thick lines) at two 
different times from the numerical simulation of the expansion 
of a two-species plasma using a Boltzmann-Vlasov-Poisson 
model. Black dashed lines are analytical profiles based on 
self-similar solutions. The spectrum of light ions shows a peak 
typical of a multispecies expansion. Reprinted figure from 
[Tikhonchuk et al.\ ( [2005| ), Plasma Phys. Contr. Fusion 47, 
B869. By permission from Institute of Physics Publishing 
(2013). 



According to the above model the peak energy of L 
ions is 



El ^ ZlTh In (4V2^N/e 



(26) 



As an important indication from this model, the mass 
ratio and the relative concentration of the two species 
might be engineered to optimize the L ion spectrum. Sev- 
eral simulation studieJ^ have been devoted to this issue 
and to the modeling of observations of multispecies spec- 
tra in both planar and spherical (droplet) targets (see 
SeclTlLEl. 



4. Numerical simulations 

Already in their simplest formulation TNSA models 
are highly nonlinear and the set of available analytical 
solutions is limited. A numerical approach can be used 
to overcome these limitations and to address additional 
effects. 

Referring to the ID problem of plasma expansion, an 
hydrodynamic two-fluid approach may be used to take 
charge separation effects into account as reported by 



Mora (2003). The hydrodynamic model, however, cannot 



take into account kinetic effects such as non-Maxwell dis- 
tribution and breakdown of equilibrium conditions. To 



[ vich et aT] (|1979|l;rKovalev et a;.[ | |2002| ; [True et al.\\198l\ ; and [ and Ruhl[ (20051; Psikal et al ('2608); R obinson et aL\ (2006| l; 



See e.g. [Brady and Arber[ l|201l|; [Bran tov et aZ.||2006|| : 



[Wickens ef aL| | |i978^ [Robinsonand Gibbon, 1 ,2007^ ; and , Robinson et aL| | |2009a l 
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FIG. 18 Electric field profiles at difi'erent times from the nu- 
merical simulation of the coUisionless expansion of a slab of 
warm plasma. Thick and dashed lines show results from a PIC 
code ( |Betti et al.\ [2005 ) and a hydrodynamics code ( |Mora[ 
|2003'), respectively. Both simulations assume a 40 thick 
proton plasma slab with initial density no — 3x 10^® cm"'' and 
electron temperature Teo = 500 KeV. The inset shows the de- 
tail of the field distributions at early times, with the field in 
the PIC simulation extending over a finite distance. Reprinted 
figure with permission from 'Romagnani et al. ( |2005 1, Phys. 
Rev. Lett. 95, 195001. Copyright 2013 by the^merican 
Physical Society. 



plasma interaction and the generation of hot electrons, 
rather than imposing a priori their number and tem- 
perature. The computational challenges and limitations 
of such large-scale simulations have been discussed in 
Sec |II.D[ In addition, most PIC simulations do not in- 
clude collisions, which may play an important role in the 



transport of hot electrons through the target (Sec II.B.3 1. 
Nevertheless, PIC simulations have been vastly used as a 
valuable support in the interpretation of measurements 
of ion acceleration and were able to reproduce at least 
qualitatively several observed features of the TNSA pic- 
ture, see e.g. [Fuchs et al.\ (|2005[); |Pukhov| (|2001[); and 



Wilks et al. (2001) 



As an alternative to the PIC method, [Gibbon et al.\ 
(2004) used a gridless, electrostatic "tree" particle code 
to simulate ion acceleration from wire targets. Such a 
code has the advantages of an unlimited spatial region for 
particles and of "automatic" inclusion of collisions, at the 
cost of being purely electrostatic so that the laser-plasma 
interaction may be modeled only phenomenologically and 
magnetic field generation is not included. 



address these latter effects a numerical solution of the 
Vlasov equation for the distribution function of electron 
and ions in phase space is needed. To this aim the PIC 
approach (Sec |II.D" ) may be used. The drawback is the 
much larger computational request with respect to hydro- 
dynamics simulations. The reason is that to obtain full 
numerical convergence and accurate, low-noise results a 
very large number of particles should be used to resolve 
the strong density variations in the plasma expansion. 

In its simplest formulation the ID simulation of coUi- 
sionless plasma expansion takes a single ion species into 
account and a limited set of parameters, such as the ini- 
tial electron temperature and the initial thickness of the 
plasma; this corresponds to fix the total energy of the 
system. Such simplified simulations already reproduce 
qualitative features observed in the experiment and may 
match measured quantities such as the ion front velocity 
with a proper choice of initial parameters. As an exam- 
ple Fig|18| shows simulation results performed to support 



experimental observations by Romagnani et al. ( 2005 1 , 
using both an hydrodynamics and a PIC code. The two 
approaches use different initial conditions, i.e. a Boltz- 
mann equilibrium for fixed ions and a zero charge density 
distribution, respectively. The latter condition enables to 
resolve in the PIC calculation the propagation of the elec- 
tron front, resulting in the electric field vanishing at the 
front position and showing a strong temporal maximum 
at the earliest instants, in agreement with experimental 
observations. 

The use of supercomputers allows to perform multi- 
dimensional PIC simulations and to simulate the laser- 



D. Comparison between models and experiments 

TNSA has been deeply investigated in a very large 
number of experiments, performed, in the past decade, 
in many laser facilities all over the world. The maximum 
observed value of the ion energy f max has been probably 
the most characterizing parameter of such experiments. 
Another important feature, mostly in the light of poten- 
tial applications, is represented by the form of the energy 
spectrum. 

All this effort resulted in an extensive collection of ex- 
perimental data, against which the predicting capability 
of the TNSA theoretical models can be tested. Moreover, 
a new generation of laser facilities will be soon available, 
and it will be possible to investigate a wider range of 
experimental parameters. Therefore the challenge of sat- 
isfactorily predicting the result of a TNSA experiment, 
providing sufficiently reliable scaling laws to extrapolate 
guidelines for the future experiments, is even more im- 
portant. 

Experimentally, great effort has been put in properly 
addressing the correlation among the above mentioned 
ion properties and the main laser and target parameters. 
Due to the importance of the laser irradiance in estab- 
lishing the regime of interaction (see Secjll|, in the lit- 
erature it has become common to report the maximum 
ion (mainly proton) energy fmax as a function of this 
parameter (see e.g. Figj4]). Collections of experimental 
data for f max have been reported in several paperj^ On 



See e.g. [Borghesi et al.\ | |2008| |2006| l; [Fuchs et al] | [2006b| ; 
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the other hand, it is apparent that the irradiance is not 
the only laser parameter playing a role in determining 
^^max- In particular, it has been established by many ex- 
perimental data that, for fixed irradiance, more energetic 
pulses lead to higher fmax- Moreover, already from the 
early experiments the strong influence both of the laser 
prepulse level and duration and of the target properties 
on fmax has been evident. We will devote Sec |III.E| to a 
dedicated discussion on these topics, while here we will 
mainly point out the general difficulties which can be en- 
countered in the attempt of providing predictions of fmax 
for comparison with experimental data. 

For all the models introduced in Sec III.C f max can 
be evaluated once the required parameters of the phys- 
ical system are known. This is actually a very delicate 
and often controversial issue because the input parame- 
ters are different in number, nature and reliability. Some 
models use laser and target parameters characterizing the 
experimental set-up, which are then known or controlled 
with well-defined precision, like e.g. the mean irradiance, 
spot radius, energy and duration of the laser pulse, and 
the thickness, chemical composition and impurity proton 
surface density of the target. Other models use as param- 
eters physical quantities determined by interaction and 
transport processes, in particular hot electron properties 



(see Sec II. B. 2 ) such as conversion efficiency rjfi, tempera- 
ture T/i, density rih, cut-off energy and beam divergence 
angle ^div These quantities may in principle be mod- 
eled and/or measured but most times are not precisely 
known. Finally, some models include purely phenomeno- 
logical parameters such as the ion acceleration time tacc 



(Sec. III.C. 2 1 and numerical parameters determined by 
fitting on some set of experimental data or numerical sim- 
ulations. The experimental and theoretical uncertainities 
and the basically different nature between model param- 
eters inevitably impose some limitations to the conclu- 
sions that one could draw on the basis of a quantitative 
comparison. 

We will briefly touch this problem following the work 



by Perego et al. (2011). The descriptions which have 



been selected there are the fluid expansion models pro- 



posed by Mora ( 2003 2005 ) , the quasi-static approaches 
of Schreiber et al. |2006') and'Passoni and Lontano (2008| 
and the "hybrid" descriptions published by Albright etaL\ 



(2006) and Robinson et al. (2006) (see Sec. III.C). The 



calculations used to implement these models and eval- 
uate fmax are also summarized. To perform the com- 
parison a database containing an extensive collection of 
published experimental parameters and results, referring 
to a wide range of laser and target parameters, has been 
considered. This analysis shows that, despite all the un- 
certainties, the predictions of the TNSA models can be 



|Krushelnick et aZ.] | |2005| ; [Perego et aZ.| | |201l| l; and |Robsoii et al.\ 
< |2007[ l. 



considered quite good, and in some cases remarkable, for 
a wide range of experimental parameters. In particular, 
quasi-static models, especially the one proposed by |Pas-| 
soni and Lontano (2008), are more suitable for the pre- 



diction of fmax- These conclusions are strongly affected 
by the estimates of the required parameters, and a more 
realistic approach to evaluate these quantities could im- 
prove the predicting capability of both expansion and 
hybrid models. 



The possibility to compare model predictions with ex- 
perimental parametric studies with well defined and con- 
trolled laser conditions, aiming at providing reliable and 
clear scaling laws, can enhance significantly the effective- 
ness of the analysis. Fig{T9] shows results from a para- 
metric study of the dependence of f max on laser power 
and duration ( Zeil a/. 20101. Several other parametric 
investigations of the dependence of f max on laser pulse 
irradiance, duration, energy and fluence have been re- 



ported ( 


Flacco et al. 




2010 


Flippo et al. 


2008b 


Fuchs 


et al. 


2006b 




Nayuki et al. 


2006 


Robson et al. 


2007) 



as well as attempts in interpreting part of these findings 



2009 


2010a 


Zani et al. 


2011 
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FIG. 19 (Color online) Experimental scaling of proton energy 
cut-off with laser power and pulse duration. Red square points 
are data from experiments performed with the DRACO laser 
at FZD (Dresden), showing a linear scaling with power in the 
short pulse (30 fs) regime. Other points are data from other 
laboratories, see |Zeil et aZ.| p010[ ) for references and details. 
The fitting Hnes correspond to the static model by Schreiber] 
et al. (2006j) with different colors (labels) corresponding to dif- 
ferent values of the pulse duration ri as given in the legenda. 
Reprinted figure from |Zeil et al.\ ( |2010[ ), New J. Phys. 12, 
045015. By permission from Institute of Physics Publishing 
(2013). 
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E. Experimental optimization 

After the first years of research, the combined vigor- 
ous development both in laser technology and advanced 
target manufacturing allowed the investigation of TNSA 
exploring a continuously increasing range of laser and 
target parameters. Actually, in most cases the two sets 
are intimately related. For example, the use of "ex- 
treme" geometrical target properties, like thicknesses 
in the sub-micrometric range, requires the availability 
of extraordinarily "clean" , prepulse-free pulses to avoid 
early target evaporation and deformation. Such pulses 
can be obtained with recently developed techniques, like 



the double plasma mirror ( 


Dromey et al. 


2004 Fuchs 


et al. 


2006a 


Thaury et al. 


2007| and references therein) , 



Optical Parametric Amplification (Shah et al. 2009) 



or Crossed Polarized Wave (XPW) generation ( [Jullien 



et al. 2005 Zaouter et al.\ 2011[ and references therein) . 



1. Energy cut-off enhancement 



Mackinnon et al. { 2002 ) studied the dependence of the 
ion acceleration on the target thickness, with the aim of 
addressing the role played by the electron temporal dy- 
namics and its effect on the formation of the sheath accel- 
erating electric field. The experimental results showed an 
increase in the peak proton energy from 6.5 to 24 MeV 
when the thickness of the Al foil target was decreased 
from 100 to 3 fim. These data clearly indicate that an 
increase in the target thickness imply a lower mean den- 
sity of the hot electrons at the surface and a consequent 
lowering of the peak proton energy. 

The influence of the laser prepulse due to amplified 
spontaneous emission (ASE) on the acceleration of pro- 
tons in thin-foil experiments has been investigated in de- 



tail by Kaluza et al. (2004). In this experiment Al foils 



of different thickness (from 0.75 to 86 /im) were used 
in connection with the possibility of changing the dura- 
tion of the ASE prepulse. The results characterized an 
optimal value for the target thickness, strongly depend- 
ing on the prepulse duration, at which the TNSA pro- 
cess leads to the highest proton energies. For the thin- 
ner targets, a prepulse-induced plasma formation at the 
rear side was able to effectively suppress TNSA, in agree- 
ment with the considerations developed in Sees |II.C] and 
III. A Related experimental work, where a wide range of 
laser parameters and different target materials have been 



considered, can be found in the literature (Fuchs et al. 



2006b Spencer ei aL 2003), 



Effective suppression of the laser prepulse level, that is 
the adoption of ultrahigh laser contrast can significantly 
alter the physical picture, since ultrathin targets, down 
to the sub-/im level, can keep their integrity until the 
interaction with the main pulse. In these conditions a 
more effective acceleration process can be expected be- 
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FIG. 20 Maximum detectable proton energy as a function 
of target thickness for high-contrast (HC) and low-contrast 
(LC) conditions. Data are shown for both backward (BWD) 
and forward (FWD) directed ions, respectively, showing the 
symmetrical behavior of TNSA for HC and ultrathin targets. 
The LC results show the existence of an "optimal" thickness 
determined by the laser prepulse causing early target distrup- 
tion, similarly to [Kaluza et a/.| ( |2004| . The laser pulse had 
65 fs duration, (0.5 ^ 1) x 10"^ W cm~^ intensity, 45° inci- 
dence and P-polarization. Reprinted figure with permission 
fromlCeccotti ei aL| (|2007|), Phys. Rev. Lett. 99,185002. 



Copyright 2013 by the American Physical Society. 



cause the refluxing and concentration of hot electrons 
in a smaller volume may lead to the establishment of 
a stronger electric field and, consequently, to higher ion 
energies. These ideas have been successfully tested by 



Neely et al. (20061, where Al target with thicknesses as 
low as 20 nm have been used in combination with 33 fs 
pulses having ASE intensity contrast reaching 10"'^''. A 
significant increase of both maximum proton energy and 
laser-to-proton energy conversion efficiency was found at 
an optimum thickness of 100 nm. Similar results have 



been obtained by Antici et al. ( 2007 ) and Ceccotti et al. 



( 2007 ) . As a further interesting feature of this latter ex- 



periment, a symmetrical TNSA on both front and rear 



sides has been demonstrated, as shown in Fig 20 when a 
sufficiently high (> 10^°) laser contrast is used. This re- 
sult confirms the universality of the TNSA process, which 
may occur also at the front side (accelerating ions in the 
backward direction) if the density profile is sharp enough. 
Very recently, using a laser pulse with similar contrast, 
40 fs duration, lO^'^ W cm~^ and irradiating targets of 



800 nm thickness, Ogura et al. (2012) reported proton 



energies up to 40 MeV, the highest value reported so far 
for pulse energies below 10 J. 

Another possible strategy to exploit the effectiveness 
in the formation of the accelerating field in mass-limited 
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targets is to reduce the lateral dimensions. Numerical in- namely flat-top hollow microcone targets (Flippo et al. 



vestigations (Psikal et al. 20081 have actually shown that 2008a), which are a modification of conical targets used 



a reduced surface leads to higher densities of hot electrons 
at the rear side of target and, thus, to higher accelerating 
electric fields. Bufi^echoux et al. (2010) experimentally 



in Fast Ignition experiments (Sec V.C). The laser pulse 



confirmed these findings showing that in targets having 
limited transverse extent, down to tens of fiia, the laser- 
generated hot electrons moving with a component of the 
velocity along the lateral direction can be reflected from 
the target edges during time scales of the same order of 
the acceleration of the most energetic ions. This trans- 
verse refluxing can result in a hotter, denser and more 
homogeneous electron sheath at the target-vacuum in- 
terface. A significant increase in the maximum proton 
energy (up to threefold) , as well as increased laser-to-ion 
conversion efficiency (up to a factor 30), can be obtained 
in these conditions, as shown in Fig|2l] Similar results, 
obtained with different laser and target parameters, have 



is focused inside the target, and starts interacting with 
the walls of the cone that it grazes while focusing down 
towards the flat top section. The remarkable reported 
result, obtained with 80 J of laser energy on the Trident 
laser at LANL, is attributed to an efficient mechanism of 
electron acceleration taking place on the inner cone walls, 
named "direct laser-light-pressure acceleration" . The re- 
sulting increase in number of high energy electrons results 
in the increase of the maximum proton energy. 

The use of targets with various structure has also been 
investigated with the particular aim to increase the ion 
energy already at relatively low laser intensitie s (below 
IQi^ W cm~^), using e.g. double layer targets (Badziak 



been found by Tresca et al. (2011 ), who also measured an 



et al. 2001 ) and more recently nanowire-covered targets 



increase in the maximum energy of protons accelerated 
from the edges of the target with decreasing target area. 

Several other attempts have been made to increase ei- 
ther the energy density of the hot electrons in the sheath 
and, consequently, the maximum proton energy. Follow- 
ing from the indications of Kaluza et al. (2004), McKenna 



et al. (2008) have investigated whether there exists an 



optimum density profile at the front of the target which 
maximizes the laser absorption. The proton cutoff energy 
was increased by 25% with respect to a sharp interface 
case at "intermediate" scale length (tens of /xms). In 
such conditions, the higher conversion efficiency into fast 
electrons was attributed to the self-focusing of the driver 
pulse. Other studies of controlled prepulse effects on ion 



(Zigler et al., 20111 for which surprisingly high energies 
up to 5.5-7.5 MeV for a 5 x 10^^ W cm'^ 40 fs laser 
pulse has been reported. 



2. Source spectrum manipulation 

Various approaches have been proposed in order to ma- 
nipulate the spectrum of TNSA protons and ions, in most 
cases with the intent to obtain narrow band peaks but 
also with the aim to enhance proton numbers through- 
out the whole spectrum or in some spectral bands, as 
required by specific applications. We will firstly review 
a number of approaches in which the proton spectrum is 



acceleration have been reported by Batani et al. (2010) 



and Flacco et al. (2008). 



modified at the source, leaving to Sec |III.E!3| approaches 
which act on the proton beam post-acceleration. 

Spectral peaks can appear as a consequence of multi- 
species plasma expansion (Sec III. C. 3 1. This effect has 



Recently, an energy cut-off increase up to 67.5 MeV, 
35% higher than for flat foil shots, has been demonstrated 



been invoked to explain observations in proton beams 
from thin foils, where the peaks appear as modulation 



by IGaillard et al.\ ^0U\ using specially devised targets, of a continuum exponential spectrum ( [Alien et al.\ |2003D 




and in experiments employing droplets of heavy water, 
where peaks are observed in the deuterium spectrum 
( Ter-Avetisyan et al. 2006). Spectral peaks have been 



observed in experiments employing high-Z metallic tar- 
gets where a plastic layer (0.5 /um PMMA) was coated 



as a dot on the rear surface of a 5 /xm Ti foil (Pfoten 



Surface area (|j 



Surface area {iim') 



hauer et al.\ |2008| [Schwoerer et al.\ |2006[ ) . These results, 
obtained on the 10 TW JETI laser in Jena, were ex- 
plained on the basis of the proton depletion approach 



FIG. 21 (Color online) Experimentally observed (a) cut-off 
proton energies and (b) conversion efficiency (for > 1.5 MeV 
protons) for 2 fim thick Au targets as a function of surface 
area, evidencing the effect of electron refluxing. The laser 
pulse had 400 fs duration, 2 x 10^^ W cm~^ intensity, 45° in- 
cidence and P-polarization. Reprinted figure with permission 
from |Buffechoux et al.\^OTo\ , Phys. Rev. Lett. 105, 015005. 
Copyright 2013 by the American Physical Society. 



first suggested by Esirkepov et al. ( 2002 1 . Robinson and 
Gibbon ( 2007 1 suggested instead that the proton density 



in the multispecies plastic layer is the important factor 
in determining the appearance of the spectral peak. Ex- 
perimental implementation required the removal of the 
native contaminant layer present at the surface and re- 
sulted in peaks in the proton spectra at 2 MeV, with 
~ 10% spread and good reproducibility (Pfotenhauer 



et al. 2008 ) 
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FIG. 22 (Color online) Ion spectra from pre-heated Pd sub- 
strate ta rgets from which hyd rogen contaminants have been 



removed ( Hegelich e£ aZ. 20061. Black curve: spectrum of C^^ 



ions. Blue curve: spectrum of the dominant substrate charge 
state Pd^'^^. Green and red curves: simulated C^^ and Pd^^^ 
spectra. Grey curve: spectrum of dominant C^"'' ions from a 
heated W target. Magenta curve: C^"^ signal from a cold 
Pd target. In the cases of black and blue curves, an ultrathin 
layer of Graphite is present on the target surface, and a quasi- 
monoenergetic spectrum appears. In the last two cases (grey 
and magenta curves) the targets have a thick layer of carbon 
contaminants and do not form a monolayer source, resulting 
in exponential-like spectra. Reprinted figu re by perm ission 
from Macmillan Publishers Ltd: (Hegelich et al. 20061, Na- 
ture 439, 441. Copyright 2006. 



Another experiment also relied on the (partial) removal 
of Hydrogen contaminants from the surface of a highZ 



Palladium target (Hegelich tt al. 2006) so that protons 



did not appear in the spectrum. Instead, monoenergetic 
features appeared in the C spectrum (specifically 0*^+), 
suggesting that an ultrathin layer of Graphite is formed 
as a result of phase changes of the Carbon compounds 
in the contaminant and that all C^'*' ions from the layer 
experience approximately the same accelerating field, as 



theoretically predicted ( Albright et al. 2006 Esirkepov 



et al. 2002 1 . Fig 22 shows spectra for targets with and 
without contaminant removal, together with hybrid sim- 
ulation predictions. 



3. Beam post-acceleration and control 

Staged acceleration employing two laser pulses on two 
separate targets has also been investigated as a possible 
route to spectral manipulation of laser-driven ion beams. 
This idea relies on accelerating a TNSA beam from a 
first target, and direct it through a second foil, which 
is irradiated by a second laser pulse at the time that a 
particular group of TNSA protons crosses the foil. These 



protons should thus experience an accelerating field as 
they transit through the rear surface of the second foil 
and gain additional energy. An experiment by [Pfoten- 



hauer et al. (|2010|, also carried out on the JETI laser. 
Peaks and dips in the spectrum 
1 MeV which correlated 



has tested this idea, 
were observed at energies of 
well with the time of flight of protons reaching the sec- 
ond target as it is irradiated, showing that the field on 
the second target slightly boosts protons in a given en- 
ergy range resulting in the spectral modification. |Burza| 



et al. (2011) reported a two stages approach employing 
spherical shell targets, irradiated by a single laser pulse, 
in which protons accelerated by TNSA at the front of 
the shell experience a second accelerating field while they 
transit through the opposite side of the shell, which mod- 
ifies the high energy end of the spectrum. The field is due 
to a hot electron charge wave spreading along the target 
surface from the interaction point, as reported in several 



experiments (. 


VIcKenna et al. 


Toncian et al. 


2006 


). 



2009a 



A different type of two stage approach was tested by 



Mar key et al.\ ( |2010[ ) on the VULCAN laser. Two pulses 
of sub-ps duration were sequentially focused, with con- 
trollable delay, on the same target in order to modify 
the temporal history of the hot electron source driving 
the TNSA, as suggested originally by [Robinson et al.\ 
(2007). An optimal delay was seen to result in an in- 



crease of energy and conversion efficiency and, addition- 
ally, a modification of the slope of the spectral profile. In 
this case, besides an optimization of hot electron produc- 
tion by the main pulse in a front-surface plasma gradi- 



ent, similar to Kaluza et al. (2004) and McKenna et al. 



(2008), the authors suggest that an additional modifi- 
cation of the proton spectrum arises from the fact that 
proton acceleration by the main pulse takes place in an al- 
ready expanding multispecies, plasma sheath at the tar- 
get rear surface. Under these conditions, the electrostatic 
field peaks at the front separating protons from heavier 
species, and re-accelerates mainly the lower energy part 



of the spectrum. Similarly, in a recent experiment, Dollar 



et al. (2011 ) have obtained spectral modifications, result- 
ing in the appearance of narrow band spectral peaks at 
^2 — 3 MeV energies, by focusing a prepulse (10~^ of 
the 10^^ W cm~^ peak intensity) on ultrathin foils a few 
tens of ps before the peak of the main pulse. 

A staged technique which acts on the protons post- 
acceleration, but employing all-optical means has been 
demonstrated by Toncian et al. (2011 2006). A tran- 



sient electric field is excited at the inner surface of a 
metal cylinder (having ~ mm diameter and length) ir- 
radiated on the outer surface by a high-intensity laser 
pulse while a laser-driven proton beam transits through 



it, see Fig 23 A-C). The field acts on the protons by mod- 
ifying their divergence leading to a narrow, coUimated 
beamlet. As the field is transient, typically lasting for 
~ 10 ps, it affects only the component transiting through 
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FIG. 23 (Color online) A-C) Schematic of laser-driven electrostatic lens. D-E) RCF stack beam profiles for protons of 9 and 
7.5 MeV, respectively, showing that the 7.5 MeV protons are focused by the fields inside the cylinder and form a black spot 
on the RCF. F): Proton spectra. Green line: spectrum obtained under same triggering conditions as in E). Black line; typical 
exponential spectrum obtained when cylinder is not triggered. Reprinted figure from Toncian et al. (20061, Science 312, 410. 
Reprinted with permission from AAAS, 2013. 



the cylinder within this time window, affecting only pro- 
tons within a narrow energy band and leading to a spike 



in the energy spectrum, as clearly visible in Fig 23 F), 
showing a 0.2 MeV band at ~ 6 MeV. Further experi- 
ments have shown that the position of the spectral peak 
can be controlled by varying the delay between the two 



laser pulses (Toncian et al. 2011 ) and confirmed that the 



focusing is chromatic, i.e. the focal position varies with 
proton energy. 

A similar approach, but employing a single pulse, was 



developed by Kar et al. ( 2008b ) for reducing the proton 



beam divergence. Also conceptually similar to the above 



described approach by Burza et al. (2011), the scheme 



employs specially designed targets in which a thin foil 
is inserted in a thicker frame, so that the charge wave 
expanding outwards from the acceleration region at the 
rear of the foil generates on the frames surface an electric 
field transverse to the expanding beam, which partially 
constrains its natural divergence. 

Other proposed methods of optical control of proton 
beam properties include beam steering triggered by shock 



waves deforming locally the target surface (Lindau et al. 



2005 Lundh et al. 2007), an effect also reported by Zeil 



et al. (20101, and control of the beam homogeneity and 



cross section profile by focusing an annular beam around 
the high-intensity interaction region, which modifies the 
properties of hot electrons refluxing through the target 



(Carroll et al. 2007). 



The high degree of beam laminarity, and the fact that 
ion emission is substantially normal to the target sur- 
face, led early on to the suggestion that by appropri- 
ately shaping the surface it should be possible to focus 



down ballistically the protons to a tight spot ( Ruhl et al. 
200T| [Wilks et al.\ [200l| ), ideally recovering the proper- 
ties of the virtual source. The idea is consistent with 
(and complementary to) observations of TNSA ions from 
wire targets, where the curvature of the target leads to 



a highly diverging beam ion with the form of an ex- 
panding disk (|Beg et al.\ |2004[ |Patel et al.\ |2003[). An 



indirect experimental demonstration of focusing was ob- 
tained via enhanced, localized heating of a secondary tar- 



get, as will be discussed in Sec V.B (Patel et al. 2003 



Suavely et al. 2007). Recently, a more direct demon- 



stration of proton beam focusing has been obtained by 
mesh projection methods in experiments where, employ- 
ing thick (250 /im) targets with hemicylindrical shape 
( Kar et aL\ [2011 ), beam focusing (down to an estimated 
25 /Ltm spot) over the whole spectrum (up to 25 MeV) was 
demonstrated. The data have highlighted the achromatic 
nature of the focusing at the different energies, consistent 
with the energy dependent variations in divergence from 
a planar foil. 

Several groups have implemented conventional acceler- 
ator techniques for energy selection or transport of laser- 
accelerated protons, in view of possible downstream ap- 
plications of the proton beam (see Sec|v]). Besides simple 
energy selection with bending magnets, the range of op- 
tions explored includes the use of pairs of quadrupole 
magnets for refocusing protons at distances in the 5- 
60 cm range and in ~ 100 /Ltm spots (Nishiuchi et al. 



2009 Schollmeier et al. 2008) or to colhmate (Ter- 



Avetisyan et al. 2008) protons within a given spectral 



band, up to 14 MeV as found by Schollmeier et al. { 2008 1 . 



A crucial parameter in this approach is the acceptance 
angle of the quadrupole system, which may limit the 
number of particles that can be focused. Large accep- 
tance pulsed solenoids (~ 9 T) were also used (Harres 
et al. 2010 Roth et al. 2009| for coUimation and trans- 
port of a large number of ^ 10^^ particles. 

The use of synchronous RF fields for phase rotation re- 
sulted in the appearance of multiple peaks across a broad- 



band spectrum (Ikegami et al. 2009). Although demon- 



strated only at relatively low energy and over low energy 
bands, this technique is in principle interesting as, rather 
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than "slicing" a portion of the spectrum, which is effec- 
tively what is done in several of the methods above, it can 
concentrate in a narrow spectral band protons originally 
contained within a larger spectral region. The phase ro- 
tation can be accompanied under the right conditions by 
a collimation effect. 

Some of these techniques have already been imple- 
mented sequentially in test beamlines operating at 1 Hz 
repetition ( Nishiuchi et al. |2010a I with a view to future 
biomedical applications (Sec V.D|. 



IV. OTHER ACCELERATION MECHANISMS 
A. Radiation Pressure Acceleration 

Electromagnetic (EM) waves carry momentum, which 
may be delivered to a non-transparent (either absorbing 
or reflecting) medium. This is the origin of radiation 
pressurcp^ whose expression for a plane, monochromatic 
EM wave of intensity / and frequency w normally incident 
on the plane surface of a medium at rest, the radiation 
pressure P is given by 



= {l + R-T)- = {2R + A)- 
c c 



where i?, T and A are the reflection, transmission and 
absorption coefflcients (with i? + T -I- A = 1) as deflned, 
e.g., in the derivation of Fresnel formulas (Jackson 1998 ) 



as a function of the refractive index and thus of the wave 
frequency, i.e. R = R{uj) and so on. Radiation pressure 
is related to the total steady ponderomotive force (PF) 



on the medium (see Sees II. A and II.B.l). The PF ef 



fectively acts on the electrons being proportional to the 
inverse of the particle mass. At the surface of an over- 
dense plasma the electrons are pushed inwards by the 
PF, leaving a charge separation layer and creating an 
electrostatic, back-holding field that in turn acts on the 
ions and leads to acceleration. 

In the case of normal incidence of a plane wave on a 
flat surface the PF density is the cycle-averaged value of 
the J X B force. In the following discussion of Radia- 
tion Pressure Acceleration (RPA) we refer to such case 
unless otherwise stated and consider only the steady ac- 
tion of radiation pressure. As discussed in Sec II.B.l 



the oscillating component of the J x B drives a sweep- 
ing oscillation at 2uj of the density profile and causes 
strong absorption and hot electron generation, except in 
the case of circular polarization for which the oscillating 



The electromagnetic theory of radiation pressure is due to James 
Clerk Maxwell I Maxweli| 1873[ l. It is however interesting that 



the Italian physicist Adolf'o Bartoli also obtained independently 
Maxwell's result in 1875 from thermodynamic considerations 
| |Barton||l884| . 



component vanishes. In the latter case, on the time scale 
of ion motion it may be assumed that the electrons are 
mostly in a mechanical equilibrium so that the PF and 
electrostatic force locally balance each other. 



1. Thick targets. Hole Boring regime 

The intense radiation pressure of the laser pulse pushes 
the surface of an overdense plasma inwards, steepening 
the density profile. For a realistic laser beam of finite 
width, the radiation pressure action drives a parabolic- 
like deformation of the plasma surface allowing the laser 
pulse penetrating deeply into the target; this process is 
commonly named "hole boring" (HB), even when refer- 
ring to a planar geometry, and it is associated with ion 
acceleration at the front side of the target. Notice that 
in the literature different definitions, such as "sweep- 



ing acceleration" ( Sentoku et al. 2003 ) or "laser piston 



(Schlegel et al.\ |2009[ ) are also used to refer to essentially 
to the same process. 

The recession velocity of the plasma surface, also 
named the HB velocity Vhb, may be simply estimated 
by balancing the EM and mass momentum flows in a 



(27) planar geometry (Denavit 1992 Robinson et al. 2009b 



Schlegel et al. 2009 Wilks et al. 1992). In the in 



stantaneous frame where the surface is at rest we ob- 
serve incoming ions with density riijhb and velocity —v^b 
bouncing back at the surface. The EM momentum flow, 
i.e. the radiation pressure, must then balance a momen- 
tum flow difference equal to ni^h{,(2mijhbVhb)vhb with 



Ihb = (1 



^hb 



-1/2 



In this reference frame, the ra- 
diation pressure is P^^^ — (2//c)(l — Vhb/c)/{l + Vhb/c) 
as can be demonstrated by a Lorentz transformation^ 
The global momentum balance thus givea^ 



2/ 1 - Vhb/c 

c 1 + Vhb/c 

Solving for Vhb this yields 
Vhb ^ hV^ 

c i + ni/2' 



nijhbi2mijhbVhb)vhb- 



(28) 



H 



Z Tir m 



miriiC'^ A n„ m. 



e 2 



«o-(29) 



^"^ Notice that the "relativistic" correction is equivalent to account 
for the energy depletion of the incident radiation in the adia- 
batic approximation. This can be easily shown by the heuristic 
model of radiation pressure as resulting from the reflection of a 
number N (per unit surface) of photons with energy-momentum 
{Hlo , -khuj / c) contained in a short bunch of duration r, corre- 
sponding to an intensity / = Nhui/r. If the surface is moving 
at velocity V = /3c, the frequency of the reflected photons is 
(jjr = Ld{l — /3)/(l + /3) and the reflection time is Tr = r/(l — 
The resulting pressure is P = |Ap|/At = {Nh/c){LU + ujr)/Tr = 
(2//c)(l-/3)/(l + /3). 

For simplicity we assume / to be independent of time. General- 
ization to a time-dependent profile I(t) is discussed by |Robinson| 
I et a/.| ( |2009b| |. 
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FIG. 24 (Color online) Hole boring acceleration by a CO2 
laser pulse in a gas jet. Left frame shows ion spectra for vari- 
ous values of the intensity /15 (in units of 10^^ W cm~'^) and 
the electron density n = ne/nc. a) /15 = 6.4, n = 6.1; b) 
/i5 = 5.5, n = 6.1; c) /15 = 5.9, n = 7.6; d) /is = 5.7, 
n = 8.0. Right frame shows the observed scaling of ion ener- 
gies with the ratio 47/nc. Reprinted figure with permission 
( |2011[ ), Phys. Rev. Lett. 106, 014801. 



from Palmer et al. 



Copyright 2013 by the American Physical Society. 



The fastest ions are those bouncing back from the surface 
in the moving frame, resulting in a maximum energy per 
nucleon in the lab frame 



n 



i + 2ni/2- 



(30) 



In the non-relativistic regime where 11 <C 1 and v^b ^ c, 
we obtain Vhb/c ~ H^/^ and E^^^ ~ 2mpC^n. 

Essentially the same results are obtained by a dynam- 
ical model of ion acceleration in the charge separation 
region at the surface ( Macchi et al. 2005 1 . Such model 



and related PIC simulations show that the ions pile up 
at the end of the skin layer producing a sharp density 
spike and causing hydrodynamical breaking and collapse 
of the electron equilibrium. This process leads to the 
production of a narrow bunch of fast ions at the velocity 
2vhb that penetrates into the plasma bulk. Eventually 
the quasi-equilibrium condition is established again and 
the process repeats itself until the laser pulse is on. HB 
acceleration is thus of pulsed nature, although on the av- 
erage it ma y be described by a steady model ( Schlegel 
[2009[ )p] 



et al. 



Eq.( 30 ) indicates that with present-day intensities high 



energies may be obtained via HB acceleration if the den- 
sity can be reduced to be slightly above ric, which is pos- 
sible if a gas jet target and a long wavelength laser, i.e. 
CO2, is used. This scheme would be interesting for appli- 
cations since it allows control of the background density, 
use of a pure proton target and high repetition rate since 
the gas is flowing. In a recent experiment [Palmer et al.\ 



^® For theoretical or simulation studies of HB by circularly polarized 
laser pulses se e also, e.g., [Chen et al. (|2008b[|; [Liseikina and] 



I Macchil ||2007|l; |Liseykina et aL| (|2008[l; [ Naumova et aL| ||2009lr 



and 



et al. 



Yin et al.' ( 2008) for single ion species case, and |Robinson| 



|| |2009c^ and |Zhang et gf., 1 ,2009^ for two ion species plasmas. 



(2011 ) employing 10 ^m wavelength, ~ 6 x 10"'^^ W cm~^ 

0.5) and a hydrogen gas 
protons of energy up 



circularly polarized pulses (oq 
jet with density of a few times ric 
to 1.2 MeV and a narrow energy spread were observed 



(Fig 24 1 . The observed ion energies were fairly consistent 
with a linear scaling with I /n^ as predicted by the HB 
model. The energies were actually higher than expected 
taking the vacuum laser intensity, suggesting that self- 
focusing in the underdense region could have increased 
the intensity in the plasma. 

We notice that Palmer et al. (2011) reported on "pro- 



tons accelerated by a radiation pressure driven shock", 
similarly to several authors who refer to HB or "piston" 
acceleration in thick targets as acceleration in the electro- 
static shock sustained by the laser pressure at the front 



surface (Schlegel ei aL 2009 Zhang ei aZ. 2009 2007b I. 



In the context of ion acceleration by laser, we prefer 
to reserve the term "shock" for the regime described in 



Sec IV. B which implies the generation of a "true" elec- 
trostatic shock wave, able to propagate into the plasma 
bulk and drive a ion acceleration there. From the point 
of view of fluid theory, a shock wave launched with some 
velocity Uaho requires the sound speed, and thus the elec- 
tron temperature, to be hot enough to prevent the Mach 
number M 

— ^sho/^s from exceeding the critical value 
M„ ~ 6.5 above which one does not have a shock but 
a "pure piston" (Forslund and Freidberg 1971). Thus, 



formation of a "true" , high speed shock wave may be in- 
hibited for circular polarization because of the reduced 
electron heating. 

Experimental evidence of HB acceleration in solid tar- 
gets is less clear at present. iBadziak et al. ( 2004 ) reported 



a series of observations of high-density, ~ keV energy ion 
pulses (plasma "blocks") for sub-relativistic irradiation 
(< 10^^ W cm~^) of sohd targets (but in the presence 
of significant preplasma). These results were interpreted 
using a model of "ponderomotive skin-layer acceleration" 
at the critical surface, a concept that sounds rather sim- 
ilar to HB-RPA. lAkH et al.\ (|2008|) reported on heating 



of solid density matter due to laser-driven density pro- 
file sweeping and shoc k formation at inte nsities up to 
lO^o W cm" 



and 



Henig et al. (2009a) reported on 



t) X 

ion acceleration by a converging shock in spherical targets 
irradiated at 1 x 10^° W cm~^. For both these latter ex- 
periments, the analysis of data and supporting PIC seems 
also compatible with HB-RPA occurring at the front sur- 
face, although the electron heating due to the use of lin- 
ear polarization complicates the picture. Indications of 
strong radiation pressure effects were also obtained from 
the modeling of coUimated, high-density plasma jets at 
the rear side of targets with a few microns thickness , at 
intensities up to 3 x 10^^ W cm~^ (Kar et al. 



2008a). It 



may be noticed that although the scaling of Eq.( 30 ) leads 



to relatively modest energies in solid-density targets, the 
foreseeable values are of interest for applications requir- 
ing large number of ions at energies of only a few MeV 
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(see SecM 



The 7 factor as a function of time is given by 



2. Thin targets. Light Sail regime 

Hole Boring RPA applies to a "thick" target, i.e. much 
thicker than the skin layer in which ion acceleration by 
the space-charge field occurs. The laser pulse penetrates 
deeper as far as adjacent surface layers are pushed into 
the target by a repeated cycle of ion bunch acceleration. 
The situation changes when a target is thin enough that 
all the ions are accelerated before the end of the laser 
pulse, i.e. a complete hole boring occurs. In such a case, 
the laser pulse is able to further accelerate ions to higher 
energies since the ions are not screened by a background 
plasma anymore. 

The thin target regime of RPA has been named "Light 
Sail" (LS) as the term is appropriate to refer to a thin 
object of finite inertia, having large surface and low mass, 
so that it can receive a significant boost from radiation 
pressure. The invention of the laser soon stimulated pos- 
sible applications of the LS concept, including visionary 



ones such as laser-driven spacecraft propulsion ( Forward 



1984). To support this idea Marx (1966) used calcula- 



tions based on the simple model of a flat, perfect mirror 
boosted by a plane wave. The analytical solution and 



scaling laws provided by such basic model ( Simmons and 



Mclnnes 1993 ) are very useful to illustrate the most ap- 



pealing features of LS-RPA, such as high conversion effi- 
ciency in the relativistic limit and the possibility to reach 
very high energies with foreseeable laser and target tech- 
nology. 

The equation of motion for a moving target ( "sail" ) in 
the laboratory frame can be obtained with the help of a 
Lorentz transformation, similarly to Eq. ( 28 ) . Neglecting 
absorption for simplicity (A = 0) we obtain 



dt 



13 



1 + /3' 



dX 
~dt 



= /3c, (31) 



where X is the position of the sail, (3 ^ V/cis its velocity 
in units of c, 7 = (1 — P^)'^^^, a = mitiil is the mass 
density per unit surface, and w' — uj[{1 — /3)/(l -I- 13)]^^"^ 
is the EM wave (laser) frequency in the rest (sail) frame. 
Notice that the intensity / is in general a function of the 
retarded time <„t = t — X/c. 



j(t) = sinh(u)H ^— -, u = -asinh(3m-t-2), (32) 

4smh(M) 3 

where fl = [ZmeO^/ AmpC^)uj and C has bee n de fined in 
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a). 



Eq.ljH]). Asymptotically, 7(t) ~ (317i)i/3 (Fig|: 

The most significant quantities can be obtained for an 
arbitrary pulse shape lit) as a function of the dimension- 
less pulse fiuence J- (the pulse energy per unit surface): 



I{t')dt'. 



(33) 



The sail velocity /3, the corresponding energy per nucleon 
E = mpCp{j — 1) and the instantaneous efficiency 77 (i.e. 
the ratio between the mechanical energy delivered to the 
sail and the incident pulse energy are given by 



/3(U) 



[1+-F(U)]'-1 



[1+Tit 
£(i„t) = nipc 



2 ^Ht...) 

2[J-(t..J + 



11' 



r/(U) 



2/3(t, 



1 + /3(U) 



1 



(34) 
(35) 
(36) 



Thus, 77—5-1 when /3(t,.ct) ^ 1- The final energy per 
nucleon £„^^^ is obtained from the total fluence Joo = 
^{tvct = 00). For a constant intensity Fao = where 
Tp is the duration of the laser pulse. In practical units 
Fao ~ '2,.2Fif,^p^^t[Q where i^ies is the fiuence in units 
niUi/l g cm~^ and iiQ = i/10 nm. 



of 10« J cm-^ pi = 
The scalings for £,„^^ are summarized in Fig 25 b). With 
present-day or near-term laser technology, fiuence values 
of 10^ J cm~'^ seem affordable, while target manufactur- 
ing can produce films of few nm thickness, e.g. Diamond- 
Like Carbon (DLC) foils. These values yield J-ao > 1 al- 
lowing to approach a regime of high efficiency, relativistic 
ions, and favorable scaling with the pulse energy. 

The above estimates have been obtained assuming a 
perfectly reflecting sail {R = 1) that, for a given surface 
density parameter C limits the laser amplitude to gq < 
C due to the onset of relativistic transparency [Eq.([5|] 
that reduces the boost on the foil. This effect suggests 
gp = C as an "optimal" condition for LS acceleration 
Tripathi et~aL\ [2009'Q that might 



(Macchi et al. 



2009c 



be however relaxed by the effect of frequency decrease in 
Analytical solutions to Eqs.pTl) exist depending on the moving fofl frame, increasing i?(a;') [see Eq.( 31 )]. For 



suitable expressions for R{uj), the simplest case being 
that of a perfectly reflecting mirror (i? = 1) and a pulse 
of constant intensity / ( Simmons and Mclnnes 1993 1 ^ 



flo > Ci all electrons are pushed away from the foil. In this 



For a constant intensity /, Eqs.l |3l[ l are identical to those for a 
charge accelerating during the Thomson scattering from a plane 
wave: see |Landau and Lifshitz] | |1962[ l who leave the solution as 
an exercise for the reader. 



■^^ The expression for 77 also follows from "photon number" con- 
servation and frequency downshift (see Section |lV. A. 1| . In the 
reflection of N photons from the mirror, the energy transferred 
to the mirror is Nn{i^ - u)r) = [2/3/(1 -|- l3)]Nhuj = T]{Nhw). 

■^^ Some authors give a similar condition for the optimal thickness 
but with slightly different numerical factors jji et aL\ |2009| |Yan| 

I et aZ.||200"8||2009a| ) 
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regime the ions in the foil undergo a Coulomb explosion 
producing a broad ion spectrum. In a composite target 
the ion field after electron expulsion might be used for 



et al. 


2011 


Qiao et al. 


2010 


2009 


Yu et al. 


2010b 



Several multi-dimensional simulation studies suggested 
to use flat-top transverse profiles to keep a quasi-lD ge- 



et al. 2008 Grech et al. 2009 ) 



monoenergetic acceleration of a proton layer (Bulanov ometry (Klimo et al. 2008 Liseykina et al. 2008 Qiao 



The interest in the LS regime was greatly stimulated 
by three-dimensional PIC simulations of thin foil accel- 



eration by Esirkepov et al. ( 2004 ) which showed that the 



temporal dependence and typical values of the ion energy 
were well described by the LS model. The simulations 
assumed a laser pulse with peak amplitude Qq = 316 
(JA^ = 1.4 X 10^^ W cm~^) and 8 cycles duration, and a 
proton slab of density 49nc and lA thickness. Most of the 
ions in a thin foil target are accelerated coherently up to 
relativistic energies (~ 1.5 GeV) as shown in Fig 26 Ac- 
cording to Esirkepov et al.\ (2004), in order for RPA to 
become the dominant acceleration mechanism the ions 
have to acquire relativistic energies already within one 
laser cycle, so that they can promptly follow electrons 
which are displaced in the longitudinal direction by the 
ponderomotive force. Later theoretical studies of such 
so-called Radiation Pressure Dominant (RPD) regime in- 



clude Rayleigh- Taylor-like instability of the foil ( jPegoraro 



and Bulanov 2007) and the effects of radiation friction 



which play a significant role at ultrarelativistic intensi- 



ties (Tamburini et al. 2010). Of particular interest is the 



possibility of a self-regulated regime where the transverse 
expansion of the foil decreases the density along the axis 
(while the frequency downshift in the foil frame compen- 
sates the effect of decreasing <^ on R{oj')), allowing for 
an increase of the ion energy at the expense of the to- 
tal number of accelerated ions ( [Bulanov e^"ar||2010a|b| . 
For a 3D expansion, theory predicts an asymptotic scal- 
ing with time of kinetic energy K{t)/mc^ ~ (3^)3/5 that 
is more favorable than for plane acceleration. This effect 
has been recently confirmed by 3D simulations ( |Tam 



burini et al. 2012") showing a higher peak energy than 



found in lower dimensionality simulations. 

The ultra-high intensities needed for RPD acceleration 
are still above present-day laser technology. However, af- 



ter the proposal of Esirkepov et al. ( 2004 ) it was realized 



that exploring the concept using pulses with circular po- 
larization (CP) at normal laser incidence would enable an 
investigation of a RPD regime at lower intensities as the- 



gets. Three papers (Klimo et al. 2008 Robinson et al. 



oretically discussed by Macchi et al. ( 2005 ) in thick tar- 



2008 Zhang et al. 2007a) independently showed that 
the use of CP allowed an optimal coupling with an ultra- 
thin foil target as well as rather monoenergetic spectra. 
Much theoretical work has been then devoted to LS-RPA 
with CP pulses, unfolding a dynamics that is much richer 
than what is included in the simple "accelerating mirror" 
model. In particular, formation of a monoenergetic ion 



distribution is not straightforward ( Eliasson et al. 2009 



Macchi et al. 2010 2009c) and may require to control 



or engineer both the pulse and target properties ( Grech 



et al. 2009 Robinson et al. 2008 ) in order to avoid tar- 



get bending that would favor electron heating, to prevent 
early pulse breakthrough due to transverse expansion and 
to keep a monoenergetic spectrum against the inhomo- 
geneous distribution of the laser intensity; for this last 
issue, a target with modulated surface density has been 
also proposed (Chen et al. 2009b). In contrast to these 



studies Yan et al. ( 2009b ) use a Gaussian intensity pro- 
file and find the formation of a narrow, high energy ion 
bunch via a self-organization mechanism somewhat sim- 



ilar to that inferred by Bulanov et al. (2010b). Another 



open issue is the stability of the foil against transverse 
perturbation^^ Very recent simulation studies charac- 
terized regimes of efficient LS-RPA for li nearly polarized 
pulses at irradiances ~ 10^^ W cm^/xm^ (Dover and Na- 



jmudin 2012 Qiao et al. 2012) 



Possible indications of the onset of LS regime have 
been recently provided in an experiment performing us- 
ing 800 fs, 3 X 10^° W cm"^ high contrast (10^) pulses 
from the VULCAN laser and very thin (^ 0.1 fim) metal- 



lic targets (Kar et al. 2012). Narrow-band spectra with 
peak energies up to ~ 10 MeV per nucleon were observed 
both for proton and heavier ions Z/A = 1/2 ions present 
as surface impurities [Fig[27]a)], while heavier bulk ions 
had a broad spectrum at lower energies. The peak ener- 
gies scaled with the fluence parameter as ^ 



ant 



[Fig 27 b)], in agreement with Eq.(35) for non-relativistic 
ions, and differently from scalings as ag, Oq or a^tp which 
have been inferred for TNSA or for other mechanisms ef- 



fective for ultrathin targets (see Sec IV. C). The Z/A = 1 
peaks are at slightly higher energy than the Z/A — 1/2 
ones, suggesting that the LS stage is followed from a mul- 



tispecies expansion (Sec III.C.3 ) in the sheath field where 
protons gain additional energy and the spectral peak sep- 
aration may be further enforced. 



The scaling plot in Fig 27 b) also contains data from 
Henig et al. p009cj who investigated LS using 45 fs, 
CP pulses at ultrahigh contrast (~ 10^-^) and ^ 5 x 
10^^ W cm~^ intensity, and few-nm DLC foils. Experi- 
mental spectra of fully ionized C®^ ions show a difference 
between linearly and circularly polarized pulses, with a 
broad peak at ~ 30 MeV appearing in the latter case, 
and reduced electron heating for CP. More recent ex- 



perimental data by Dollar et al. (2012) using tightly fo- 
cused (//I) pulses with intensity up to 2 x 10^^ W cm^'^ 
showed weak difference between CP and LP, which was 
attributed to the early deformation of the thin targets 



23 See e.g . jAdusumilli ef"aL](|2012t;|Chen et aLl|201l||;|Klimo et al.\ 
(|2008^; IPegoraro and Bulanov| ( |2007| l; |Tikhoiichuk| | |2010[ |; aiid 
|Yu ei af.| ( |20T0b| 
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causing excessive electron heating. Very recently, pre- 
liminary indications of a transverse instability resulting 
in spatial modulations of the accelerated proton beam 



have been reported (Palmer et al. 20121 



To summarize the experimental evidence, so far there 
is a fairly clear confirmation of the expected LS scaling, 
but also indications of significant detrimental effects. The 
observed ion spectrum is relatively broad, suggesting that 
transverse inhomogeneity and heating effects need to be 
reduced. In perspective, the relatively slow scaling of 
the energy gain with time might pose the challenge to 
increase the acceleration length against the effect of, e.g., 
pulse diffraction and instabilities. 



B. Collisionless Shock Acceleration 



strongly relativistic limit ap ^ 1 the condition to obtain 
radiation-pressure driven supersonic shocks {M > 1) can 
be written as a/Soo > ne/uc- The reflected ions may get 
further acceleration by the transient sheath field at the 
rear surface as in TNSA, eventually producing a plateau 
in the ion spectrum. A similar signature was observed 



experimentally by Zepf et al. { 2003 1 and thus intepreted 



as evidence of front side contribution to ion acceleration, 
in contrast to pure TNSA at the rear side of the tar- 
get. In particular conditions, the staged CSA-TNSA ac- 
celeration might produce the highest energy component 
in the ion spectrum as observed in simulation studies 



(Chen et al. 2007 d'Humieres et al. 2005) which how- 



ever also suggest lower efficiency and brilliance with re- 
spect to pure TNSA. 



Acceleration of particles by shock waves (briefly, 
shocks) in plasmas is a problem of central interest in as- 
trophysics (Martins et al. 2009). The existence of an ion 



component that is reflected by the shock front is actually 
integral to the formation of the collisionless, electrostatic 
shock waves in basic fluid theory where the electrons are 



assumed to be in a Boltzmann equilibrium ( 


Forslund and 


FreidbergI 1971| |Forslund and Shonk| |1970t 


Tidman and 


Krall| 1971). In the frame moving at the shock veloc- 



ity, ions are reflected by the shock if the height of the 
electrostatic potential barrier ^^^x at the front is such 
that Ze<I>,-„ax > ■mivl/2, being vi the velocity of the ion 
component in the shock frame. Behind the shock front, 
the fields have an oscillatory behavior. Reflected ions ini- 
tially at rest acquire a velocity in the lab frame equal to 
2i'Bho where Wgho is the shock front velocity. 

Collisionless shock acceleration (CSA) was proposed as 
an ion acceleration mechanism in superintense laser in- 
teraction with an overdense p lasm a on the bas i s of P IC 
simulations by |Denavit| ( |1992| and |Silva erZ] ( |2004| 



In the latter work simulation showed the generation of 
shocks with high Mach numbers M = Weho/cs = 2 — 3, 
where the sound speed is estimated using for the "tem- 
perature" the hot electron energy, i.e. Th — £p ([6|. The 
shocks are generated at the front surface with a veloc- 



ity close to Vhb given by Eq.(29l, consistently with the 
assumption that they are driven by the piston action of 
radiation pressure. By estimating 



Vhbr] in the 



In experiments on underdense plasmas created either by using 
gas jet targets ( |Wei et aL\\2004l or by the effect of long prepulses 
in solid targets l |Habara et al. \ |2004b[ |, the observation of ion 
acceleration along the radial direction has been attributed to 
radial shock generation in a laser-driven channel. 
Notice that tigho — ^hb implies that "reflected" ions directed into 
the bulk will have a velocity ~ 2vf^i,, i.e. twice the surface reces- 
sion velocity, as the fastest ions generated by the piston action 
in HB acceleration (Sec |IV.A.l[ This similarity may explain why 
HB and CSA are often confused in the literature. 



Very recently, CSA has been indicated as the mecha- 
nism responsible for monoenergetic acceleration of pro- 
tons up to 22 MeV in the interaction of CO2 laser 
pulses with Hydrogen gas jets at intensitie s up to 6.5 x 
10^^ W cm~^ corresponding to oq = 2.5 (Haberberger 



et al. 2012). The particular temporal structure of the 
laser pulse, i.e. a 100 ps train of 3 ps pulses, was found 
to be essential for the acceleration mechanism, since no 
spectral peaks were observed for a smooth, not modu- 
lated pulse. Comparison with PIC simulations suggested 
that the multiple pulses lead to efficient generation of 
suprathermal electrons, and that the latter (rather than 
radiation pressure) drive the shocks which eventually ac- 
celerate protons. Simulations also suggest that the pro- 
cess could scale in order to produce 200 MeV protons 
at 10^^ W cm~^, that may foreseeable with future CO2 
laser development. Such scheme based on gas lasers and 
gas jet target would have the remarkable advantage of 
high-repetition rate operation, but the efficiency per shot 
might be low with respect to other approaches: in the ex- 



periment of Haberberger et al. (2012) the number of ions 
(~ 2.5 X 10*^ ) in the narrow spectral peak at ~ 22 MeV 
for a 60 J pulse energy implies a conversion efficiency of 
^ 10-8. 



In addition to collisionless shocks, the standard fluid 



theory also predict solitons (Tidman and Krall 
propagating at the velocity v, 



1971[ ) 



with 1 < v^ai/cs 
These solitons are characterized by Ze$„„x < ^iv'^ai/'^ 
and are thus transparent to background ions "by con- 
struction". However, generation of electrostatic solitons 
may lead to ion acceleration in some circumstances, e.g. 
when the soliton breaks in the expanding rear sheath due 



to the effect of the plasma flow ( Zhidkov et al. . 2002 1 . 



Additional simulation studies of shock and solitary wave 



acceleration are reported by He et al. (2007); Liu et al. 



(2009): and Macchi et al. (2012). 
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C. Transparency regime. Break-out afterburner 



D. Acceleration in near-critical and underdense plasmas 



If ultrathin foils are used as targets (which requires 
ultrahigh-contrast, prepulse-free conditions), the expan- 
sion of the foil may lead to the onset of transparency 
during the short pulse interaction, when the electron den- 
sity Ue is further decreased down to the cut-off value (of 



the order of jric due to relativistic effects, see Sec II. A). 
While this effect limits the energy attainable via RPA 
(Sec |IV.Al2 ), it can lead to enhanced ion acceleration 
via different mechanisms. 

Several related experiments were performed at the 
TRIDENT laser facility at Los Alamos National Labo- 
ratory (LANL), using pulse durations in the 500-700 fs 
range. C®^ ions with energies up to 15 MeV per nucleon 
were observed by irradiating DLC foils with ~ 40 J, 
~ 7 X 10^^ W cm~^ pulses, for an optimal thickness 
of 30 nm that is determined by the condition that rel- 



ativistic transparency occurs at the pulse peak (Henig 



et al. 2009b). Fig 29 shows spectra for different polar- 
izations. For more energetic and intense pulses (~ 80 J, 
10^1 W cm"^) and thicker targets (140 nm), broad 
spectra with higher cut-off energies beyond 40 MeV per 
nucleon were observed, and the inferred conversion ef- 
ficiency was ^ 10% (Hegelich et al. 2011). Narrower 
C^~^ spectra {A£i/£i ~ 15 — 20%) at lower energies 
('-^ 3 — 10 MeV) were observed using either loose focus- 



ing or circular polarization (Hegelich et al. 2011 Jung 



et al. 2011b). Very recently, energies up 80 MeV per 
nucleon for Carbon and 120 MeV per protons have been 
communicated (Hegelich 2011). The onset of relativis- 



tic transparency in these conditions has been recently 
investigated in detail with ultrafast temporal resolution 



( Palaniyappan et al. 2012) 



Simulation studies of this regime show that the in- 
crease of the cut-off energy is related to enhanced and 
volumetric heating of electrons as the target becomes 
transparent, leading to a stronger accelerating field for 
ions; the name "Break-Out Afterburner" (BOA) has been 
proposed for such regime by the Los Alamos group (Yin 



et al. 2007). 3D PIC simulations of BOA have been re- 



ported by Yin et al. (2011a I. Modeling of BOA is not 



simple as the process appears to involve different stages. 
Analytical descriptions of BOA have been reported by 



Albright et al. (2010) and Yan et al. (2010) and a scaling 



of the maximum ion energy £„ax — (1 -I- 2a)ZTe, with 
the electron temperature and a a phenomenological pa- 
rameter (estimated to be ~ 3 from simulations), has been 
proposed. A fastly growing relativistic Buneman insta- 
bility, excited due to the relative drift between electron 
and ions, has been invoked as a mechanism enhancing 
the couphng with ions (Albright et al.\ 2007). Theoret- 



ical explanations for narrow C°~*' spectra, based on an 
electromagnetic "ion-soliton" model (fundamentally dif- 



ferent from electrostatic solitons described in Sec IV. B ) 



are discussed by Yin et al. (201 lb I. 



A number of studies has been devoted to ion accelera- 
tion in "near-critical" plasmas with electron density close 
to the cut-off value (rig — rig), in order to allow a more 
efficient generation of hot electrons to drive TNSA. Pro- 
duction of a low-density plasma by a laser prepulse has 
been investigated for laser and target parameters such 
that at the time of interaction with the main short pulse 



the preplasma was either underdense (Matsukado et al. 



2003) or slightly overdense (Yogo et al. 2008); in the 
latter experiment, protons up to 3.8 MeV are observed 
at 10"'^^ W cm~^ intensity, and directed slightly off the 
normal to the target rear side. The analysis of these ex- 
periments gave indication of a regime where the pressure 
due to self-generated magnetic field at the rear surface 
strongly contributes to charge separation. 

An alternative strategy to reduce the electron den- 
sity is to use special target materials such as foams, 
which may be manufactured in order to have an av- 
erage value of rie slightly larger, or even lower than 
Ue (the average is meant over a length larger than the 
typical sub-micrometric scale of inhomogeneity) . Ex- 
perimentally, proton acceleration in low-density foams 
(rip, = 0.9 — 30ne) has been investigated by [Willingale 



et al. 



(2009 2011a) at intensities up to lO^i W cm"^. 



showing that the proton energy is close to that obtained 
for solid foils and the same laser pulse for the lowest 
density value (rig ~ 0.9nc). In this experiment, proton 
acceleration has been mostly investigated as an indica- 
tion of the onset of relativistic transparency, leading to 
enhanced laser penetration and coUimation of hot elec- 
trons and ions by self-generated magnetic fields. Recent 
simulation studies of ion acceleration in solid target cov- 



ered with foam layers have been also reported ( [Nakamura 
et al.\ |2010[ [Sgattoni et'a!] |2012| . 

Experimental investigations of ion acceleration using 
gas jet targets, with typical densities below 10^° cm~'^, 
have been also performed. These experiments include 
the already described investigations of hole boring RPA 



(Sec IV. A. 1) and shock acceleration (Sec IV. B I using 
CO2 lasers for which gas jets are near-critical targets. 
Using optical or near-infrared lasers, several experiments 
of high-intensity laser interaction with underdense gas 
jets have reported observations of energetic ions accel- 
erated in the radial direction with respect to the laser 



pulse propagation axis (|Krushelnick et al. 


1999 Sarkisov 


et al. 


1999 Wei et al. 


2004 


) by the electric field created 



by the electron displacement in the channel drilled by the 
ponderomotive force. The ion spectrum may provide in- 
formation on the self-focusing and channeling dynamics 
of the laser pulse and the acceleration mechanism shows 
indeed some similarity with those active in the interac- 



tion with sohd targets (Macchi et al. 2009a). However, 



radial acceleration of ions is of modest interest for appli- 
cations since the ions are not collimated at all. 
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A collimated emission in the forward direction from an 
underdense He gas jet was reported by [Willingale et al. 
( 2006 ) . Using a laser pulse of 1 ps duration, energy up to 



340 J and intensity up to 6 x 10^° W cm~^, He ions up to 
40 MeV were observed collimated in a beam with < 10° 
aperture. The data were interpreted by assuming that a 
large electric field was generated at the rear side of the 
gas jet by escaping hot electrons. Differently from TNSA 
in solid targets, the mechanism was considered to be ef- 
fective despite the relatively long density scalelength at 
the rear surface because a larger fraction of hot electrons 
was generated by electron acceleration in the underdense 
plasma. Further analysis of simulations of the experi- 



ment (Willingale et al. 2007) showed also a significant 



contribution due to the generation of a quasi-static mag- 
netic field at the rear surface, which in turns enhances 
the accelerating electric field via magnetic pressure and 
induction effects according to the model by |Bulanov and 



Esirkepov ( 2007 1 that was also used to explain the above 



mentioned experimental results by Yogo et al. { 2008 1 in 
a near-critical plasma. Figj30] shows a sketch of such 
mechanism. 

Magnetic-field sustained ion acceleration was also in- 
dicated as the dominant mechanism in an experiment 



by [Fukuda et al.\ ( |2009[ ), where ions in a 10-20 MeV 
range and collimated in a 3.4° aperture cone were ob- 
served in the interaction of a 7 x 10^^ W cm^^, 40 fs 
laser pulse with an underdense gas jet where CO2 clus- 
ters were formed. The role of the clusters was apparently 
that of enhancing the self- channeling and focusing of the 
laser pulse, leading to an increase of the intensity in the 
plasma, rather than contributing to ion acceleration via 
cluster explosions. 

Generation of collimated ions from underdense plas- 
mas at ultrahigh intensities (> 10^"'^ W cm~^) was in- 
vestigated theoretically and with numerical simulations 



aheady more than a decade ago (Bulanov et aL, 2000 
Esirkepov et al, 1999| Sentoku et al. 2000). In par- 
ticular in these papers it was predicted that, for oo > 
(rrii/meY^'^ — 43A^/^, the effective inertia of the highly 
relativistic electrons in the laser field becomes compara- 
ble to those of ions. As a consequence the ions closely fol- 
low the electron displacement due to the ponderomotive 
action and the acceleration process may become similar 
to what is observed in an overdense plasma. A few more 
recent simulation studies investigated a regime where a 



small ion target is placed in an underdense plasma ( Shen 



et al. 2009 Yu et al. 2010a I. The superintense laser 



pulse accelerates and overruns the ion target and then 
generates a wakefield in the underdense plasma, where 
ions may be trapped and accelerated in a way similar to 
the well-known scheme for laser acceleration of electrons 



E. Resistively enhanced acceleration 

Already during the "front vs rear side acceleration" 
debate related to experiments reported in 2000 (Secjl|, 
it was suggested that protons may also be accelerated in 
the target bulk through a mechanism depending on the 



target resistivity rj ( Davies 2002 ) . The electric field gen- 



erated in the target bulk to provide the return current 



E = jr/?7 (see Sec H.B.3 ) increases for low 77 reducing the 
penetration of hot electrons through the target and at 
the same time favoring acceleration in the front and bulk 
regions versus TNSA. The mechanism has been theoret- 



ically investigated by Gibbon (2005a) using a coUisional 
tree-code approach. 

Indications of dominant front side acceleration due to 
resistivity effects have been reported in solid plastic tar- 



gets (Lee et al. 2011 2008 1 and also in low-density foams 



( Li et al. 2005 ) where an anomalously high resistivity 



might be due to spatially localized fields in the locally 
inhomogenous material. 



V. CURRENT AND FUTURE APPLICATIONS 
A. Proton radiography 

The use of ion beams, and particularly proton beams, 
for radiographic applications was first proposed in the 
1960s (Koehler 1968). Quasi-monochromatic beams of 



ions from conventional accelerators have been used for 
detecting aereal density variations in samples via modifi- 
cations of the proton beam density cross section, caused 
by differential stopping of the ions, or by scattering. Ra- 
diography with very high energy protons 1 — 10 GeV) 
is being developed as a tool for weapon testing (King 



et al. 1999). Ion beams from accelerators have also been 



employed in some occasions for electric field measure- 
ments in plasmas, via the detection of the proton de- 
fiection, 



e.g. 



Mendel and Olsen (1975). In practice. 



(Esarey et al.. 2009). In those simulations GeV energies the target (Borghesi et al. 2004). 



the difficulties and high cost involved in coupling exter- 
nally produced particle beams of sufficiently high energy 
to laser-plasma experiments (or indeed magnetic confine- 
ment experiments) and the relatively long duration of ion 
pulses produced from conventional accelerators have lim- 
ited the application of such diagnostic techniques. 

The unique properties of protons from high intensity 
laser-matter interactions, particularly in terms of spatial 
quality and temporal duration, have opened up a totally 
new area of application of proton probing or radiogra- 
phy. As seen in Sec jllLB] the protons emitted from a 
laser-irradiated foil by TNSA can be described as emit- 
ted from a virtual, point-like source located in front of 

A point-projection 



were reached, but the required laser pulses should have 
multi-petawatt power and multi-kJ energy, that is still 
far beyond present-day laser technology. 



imaging scheme is therefore automatically achieved with 
magnification M set by the geometrical distances at play. 
Backlighting with laser-driven protons has intrinsically 
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high spatial resolution, which, for negligible scattering in 
the investigated sample, is determined by the size d of 
the virtual proton source and the width 6s of the point 
spread function of the detector (mainly due to scattering 
near the end of the proton range) , offering the possibil- 
ity of resolving details with spatial dimensions of a few 



Two main arrangements have been explored. In proton 
imaging, i.e. simple backlighting projection of the sample 



(Borghesi et al. 2002b 2001), the deflections cause local 



modulations in the proton density Up across the proton 
beam cross section, which, under the approximation of 
small deflections, can be written as 



/Ltm. As discussed in Sec II. E multilayer detector arrange- 
ments employing RCFs or CR39 layers offer the possi- 
bility of energy-resolved measurements despite the broad 
spectrum. Energy dispersion provides the technique with 
an intrinsic multi- frame capability. In fact, since the sam- 
ple to be probed is situated at a finite distance from the 
source, protons with different energies reach it at differ- 
ent times. As the detector performs spectral selection, 
each RCF layer contains, in a first approximation, infor- 
mation pertaining to a particular time, so that a movie of 
the interaction made up of discrete frames can be taken 
in a single shot. Depending on the experimental condi- 
tions, 2D proton deflection map frames spanning up to 
100 ps can be obtained. The ultimate limit of the tem- 
poral resolution is given by the duration of the proton 
burst at the source, which is of the order of the laser 
pulse duration. 

Several radiographic applications of laser-produced 
protons have been reported to date and radiographs of 
objects for various size and thickness (down to a few 
^m) have been obtained ( Borghesi et a/. [ [2004 Cob- 



eL 
' 2epM 



+b/2 



E 



-b/2 



xB]dx (37) 



where Vp and Ep are the proton velocity and energy, M 
the projection magnification, L the distance between the 
plasma and the detector, and the integral is along the tra- 
jectory of the protons, crossing a region \x\ < b/2 where 
the fields E and B to be probed are present. Under 
simplified assumptions the above formula can be used 



to yield line-averaged values of the fields (Sarri et al. 



2010a). In proton deflectometry, thin meshes are inserted 
in the beam between the proton source and the object 
as "markers" of the different parts of the proton beam 
cross sections (Mackinnon et al. 2004). The meshes im- 



press a modulation pattern in the beam before propagat- 
ing through the electric field configuration to be probed. 
The beam is in this way effectively divided in a series of 
beamlets, and their deflection can be obtained directly 
from the pattern deformation. A technique employing 
two grids to generate a set of Moire fringes has also been 
proposed as a way to increase the sensitivity to small 



ble et al.\ [2002} [Mackinnon et al.\ |2006| |Roth et aL\ electric fields ([Mackinnon TTal] [2003 1 . 



2002). Density diagnosis via proton radiography has po- 
tential application in Inertial Confinement Fusion (ICE) 
A preliminary test studying the compression of empty 
CH shells under multi-beam isotropic irradiation at the 
moderate irradiance of 10^^ W cm^^ has been studied 
in an experiment carried out at the Rutherford Apple- 
ton Laboratory (Mackinnon et al. 2006). Radiographs 



The proton probing technique has provided uniquely 
detailed information on nonlinear phenomena in high- 
intensity laser-plasma interactions, such as ion acoustic 



solitons and collisionless shock waves ( Romagnani et al. 



2008 ) , phase-space electron holes ( Sarri et al. 2010b ) , the 



charge-displacement channel formation dynamics foUow- 



of the target at various stages of compression were ob- 
tained. Modelling of proton propagation through target 
and detector carried out using Monte Carlo codes per- 
mits the retrieval of density and core size at maximum 
compression (3 g/cc, 80 /im) in good agreement with hy- 
drodynamic simulations. Radiographic analysis of cylin- 
drically compressed mat ter (|Volpe et~ar. 2011) and of 



ing relativistic self-focusing of laser pulses ( 


Kar et al. 


2007 


Sarri et al. 2010c 


Willingale et al. 


2011b 1 and the 



evolution of remnants of coherent electromagnetic struc- 



tures and instabilities of various type (Borghesi et al. 



2002a 



[20051 [RoSTagnani a/.l [20T0l [Sam et aL[|2011b 



2010d). Application to ns laser-produced plasmas of ICE 



shock wave propagation (Ravasio et al. 2010) has been 



interest has also allowed to investigate laser filamentation 



in underdense plasmas (Lancia et al. 2011 



also carried out, although the available detail with low 
energy protons was limited. 

The most successful applications to date of proton 
probing are related to the detection of electric and mag- 



et al. 



Sarri et al. 



2009 



ingale et al. 



2011a), plasma expansion inside hohlraums (Sarri et al. 
2010a) and self-generation of magnetic fields (Cecchetti 
"Nilson et a/.[ [2006j [Sarri eToIl |2011a| (Will 



2010). As an example of the use of time- 



netic fields in plasmas (Borghesi et al. 2002b Mackin- 



non et al. 2004 ) . Jointly with a parallel technique using 



resolved proton diagnostic, Eig[l3| reports data from an 
experiment where the protons are used to probe the rear 
of a foil following ultraintense irradiation of the front 



monoenergetic protons from fusion reactions driven from of the foil ( Romagnani et al. 2005 1 . The probe proton 



laser-driven compressions (Li et al. 2006), proton prob- 



ing with laser-accelerated protons has provided in this 
way novel and unique information on a broad range of 
plasma phenomena. The high temporal resolution is here 
fundamental in allowing the detection of highly transient 
fields following short pulse interaction. 



pattern is modified by the fields appearing at the target 
rear as a consequence of the interaction, and the tech- 
nique effectively allows spatially and temporally resolved 
mapping of the electrostatic fields associated to TNSA 
acceleration from the foil (see Sec III.A). Eig 13 a) shows 



the set-up for both imaging and deflectometry measure- 
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merits. Frames b)-g) of Fig 13 correspond to proton im- cal and planetary studies (Ichimaru 1982 Lee et al. 



ages at different times taken in a single shot, resolving 
the expansion of the plasma sheath and highlighting the 
multiframe capability of this diagnostic. 



2003). The high-energy flux and short temporal dura- 



tion of laser-generated proton beams are crucial param- 
eters for this class of applications. WDM states can be 



It could be noticed than on the basis of Eq.(37) it may achieved by several other means, e.g. X-ray heating ( Tal 



not be possible in principle to attribute unambigously 
the measured deflections to the sole action of either elec- 
tric or magnetic fields. Confidence in the interpretation 
of observed patterns can be increased by supporting the 
analysis method of both imaging and deflectometry data 
with particle tracing codes. Such codes simulate the 
propagation of the protons through a given space- and 
time-dependent field configuration, which can be modi- 
fied iteratively until the computational proton profile re- 
produces the experimental ones. State-of-the-art tracers 
allow realistic simulations including experimental proton 
spectrum and emission geometry, as well as detector re- 
sponse. Moreover, in some specific experiments it was 
possible to give evidence of magnetic fields, discriminat- 
ing their effect on probe protons from that due to elec- 
tric fields, thanks to the possibility of different probing 



lents et al. 2009 1 and shock compression ( Kritcher et al. 



2008). However, when studying fundamental properties 



of WDM, such as equation of state or opacity, it is desir- 
able to generate large volumes of uniformly heated ma- 
terial; ion beams, which can heat the material in depth, 
are in principle better suited to this purpose than the 
methods described above. 

Heating of solid density material with ions can be 
achieved with accelerator-based or electrical-pulsed ion 
sources, see e.g. IBailey et al.\ (|1990|); IHoffmann et al. 



(2000); and Tahir et al. (2006). However the relatively 



directions ( Cecchetti et al. 2009 ) or even exploiting the 



divergence of the probe beam (Romagnani et aL, 2010). 



An example is given in FigjST] where the presence of an 
azimuthal B-field has been revealed with mesh deflec- 
tometry by either a compression or outwards dilation of 
the mesh lines, depending on whether the B-field has 
clockwise or counter clockwise direction compared to the 
propagation direction of the probe beam. 

The divergence of the probe beam also implies that 
the effective probing time is a function of the position on 
the image plane because of the different time of fiight for 
protons at different angles. This effect has to be taken 
into account for measurement of field structures prop- 



long durations of ion pulses from these sources (1-10 ns) 
means that the materials undergo significant hydrody- 
namic expansion already during the heating period. On 
the contrary, laser-generated proton beams, emitted in 
ps bursts, provide a means of very rapid heating, on a 
timescale shorter than the hydrodynamic timescale. By 
minimizing the distance between the ion source and the 
sample to be heated, it is possible to limit the heating 
time to tens of ps. The target then stays at near-solid 
density before significant expansion occur, and the WDM 
properties can be investigated within this temporal win- 
dow. 

The first demonstration of laser-generated proton heat- 
ing was obtained by Patel et al. (2003). In this experi- 



agating at relativistic speeds (Kar et al. 2007 Quinn 



ment a 10 J pulse from the 100 fs JanUSP laser at LLNL 
was focused onto an Al foil producing a 100-200 mJ pro- 
ton beam used to heat a second Al foil Target heating 
was monitored via time-resolved rear surface emission, 
as shown in Fig|33) A focused proton beam, produced 



et al. 2009a) and actually may improve the capability from a spherically-shaped target (Sec |III.E73 ), was seen 



to characterize such structures, as it was obtained by 



a slightly modified arrangement (Quinn et al. 2009b I. 



This allowed to characterize the ultrafast, transient field 
front associated to the early stage of TNSA where elec- 



to heat a smaller region to a significantly higher tem- 
perature, 23 eV vs 4 eV) with respect to a focused 
beam. With a similar ion focusing arrangement on a 
higher energy laser system, Gekko at ILE Osaka, Suavely] 



tromagnetic effects come into play (see Fig 32 1 . A proton et al. ( 2007 ) demonstrated secondary target heating up 



streak deflectometry technique for obtaining continuous 
temporal mapping (but only one spatial dimension is re- 
solved) has also been proposed, in which the energy res- 
olution is done by means of magnet dispersion ( Sokollik 
eral||2008| . 



to 80 eV by imaging both visible and extreme-ultraviolet 
Planckian emission from the targets rear surface. 

Subsequent experiments have investigated the proper- 
ties of the WDM produced in this way with a number 
of diagnostics, either passive or in pump-probe conflg- 
urations, combined to self-consistent modeling of sam- 
ple heating and expansion. Warm solid Al at tempera- 



B. Production of Warm Dense Matter 

Laser-driven ions have found application in a num- 
ber of experiments aimed to heat up solid density mat- 
ter via isochoric heating, and create so-called Warm 
Dense Matter (WDM) states, i.e. matter at 1-10 times 



tures up to 15-20 eV (Dyer et al. 2008 Mancic et al. 



2010b) and Carbon up to ~ 2 eV (Roth et al. 



Dyer et 



2009) 



(2008) 



solid density and temperatures up to 100 eV (Koenig 



et al. 2005 ) of broad relevance to material, geophysi- 



have been produced in this manner 
reconstructed the Equation of State (EOS) of the heated 
material by measuring the temperature and expansion 
rate of the heated target, via streaked thermal emission 
and chirped pulse interferometry. 

Pump-probe arrangements have been used in recent. 
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more sophisticated experiments, which have provided to- 
tally novel information on the transition phase between 
cold solid and plasmas in isochorically heated Al and 



C targets. Mancic et al. (2010a) have investigated the 



short range disordering of warm Al at solid density by 
applying time-resolved X-ray absorption near-edge spec- 



troscopy (see Fig 34 ) . Progressive smoothing of spectral 
features near the K-edge allowed to place an upper bound 
on the onset of ion lattice disorder within the heated 
solid-density medium of 10 ps. Pelka et al. (20101 have 



recenty diagnosed ultrafast melting of carbon samples, 
by X-ray scattering techniques, which allowed determi- 
nation of the fraction of melted carbon in the heated 
sample. Comparison to predictions based on different 
theoretical descriptions of the EOS of Carbon, indicates 
a departure from existing models, with implications for 
planetary core studies. 

In all the experiments mentioned above the isochoric 
heating by the protons is volumetric, but not uniform 
(Brambrink et al. 2007), see e.g. Fig 34 c). Uniform 



heating would require some degree of proton energy se- 
lection, and choosing the sample thickness so that the 
Bragg peak of the selected protons does not fall within 
the sample, as suggested for example by |Schollmeier et af. 
(120081). 



C. Fast Ignition of fusion targets 



The traditional route to ICF (Atzeni and Meyer-ter- 



Vehn, 2004) relies on the driven implosion of a pellet of 



thermonuclear fuel (a DT mixture). Ignition occurs in a 
central "hot spot" following pulse compression. This ap- 
proach requires an extremely high symmetry and is prone 
to hydrodynamics instabilities, making ICF a historically 
difficult goal. 



In the Fast Ignition (FI) concept (see Key ( 2007 ) for a 
compact review) ignition is driven by an external trigger, 
creating the hot spot in a time much shorter than the typ- 
ical fuel disassembly time. Hence, ignition is separated 
from pulse compression. The FI approach might relax 
symmetry and stability requirements, reduce the energy 
need for ignition and allow fuel burn in a isochoric regime 
with high fusion gain. 



In the original FI proposal by Tabak et al. (1994), the 



foreseen ignitor beam was composed by multi-MeV elec- 
trons accelerated by a petawatt laser pulse via the mecha- 
nisms described in Sec |II.B[ Subsequent research showed 
that, besides generating an electron beam with enough 
power to ignite, most problematic were the issues of en- 
ergy transport and deposition in the core. Concerning 
the latter issue, the energy deposition profile of electrons 
is a smooth function, making it difficult to produce a 
localized hot spot. 

The observation of efficient generation of multi-MeV 



proton beams in Petawatt experiments (Hatchett et al. 



2000 Suavely et al. 2000 ) soon stimulated the proposal 



of the use of such protons as an ignitor beam ( |Roth et al. 
2001). The most promising features of proton beam ig- 



nition as claimed in the paper were the highly localized 
energy deposition profile (see Figjs]), the low emittance 
of the beam and its focusability, for instance by paraboli- 
cally shaping the rear side of the proton- producing target 



as suggested by numerical simulations (Ruhl et al. 2001 



Wilks et al. 2001 1 



Detailed calculations by Atzeni et al. ( 2002 ) and Tem- 



poral et al. ( 2002 1 addressed in particular the effects of 



the quasi-thermal energy distribution typical of TNSA 
protons, and of the related temporal dispersion. The 
latter could be beneficial for energy deposition since the 
proton stopping range increases with plasma tempera- 
ture. Hence, heating due to the more energetic protons 
favors energy deposition by the less energetic ones which 
arrive later in the dense fuel region. A fit of simulations 
for a proton temperature of 5 MeV provided the follow- 
ing estimate of the ignition energjj^ S^^ as a function of 
fuel density p and distance d between proton source and 
fuel core: 

~ 90(d/mm)"-^(p/100 g cm-3)-i-3 kJ. (38) 

Integration of the foil inside the cone of conically-guided 
ICF targets already designed for electron FI was then 
proposed in order to reduce d and thus fig. This raised 
the issue of shielding the foil from pre-heating caused, 
e.g., by external radiation, which may jeopardize efh- 



cient TNSA (see SecIII.E.ll; a preliminary analysis is 
mentioned by Geissel et al. ( 2005[ ). FigjSS] sketches the 
target and foil assembly and summarizes suitable param- 



eters for protons FI with cone targets. Temporal (2006) 



and Temporal et al. ( 2008 ) investigated a similar scheme 



but using two proton beams with suitably shaped radial 
profiles, obtaining a 40% reduction of the ignition energy. 



Independently from the paper by Roth et al. (20011, 



fast ignition by laser accelerated light ion beams was pro- 
posed by B ychenkov et al.^ (2001a) . The use of deuterons 
and Beryllium ions was investigated, and it was sug- 
gested that using ions with Z > 1 could be advantageous 
because of their higher stopping power with respect to 
protons. Calculations by Honrubia et al. (20091 showed 



that, differently from the proton-based scheme, ions with 
a narrow energy spread S£/£ would allow to lower the 
ignition threshold: for S£/£ = 0.1, < 10 kJ might 
be obtained. This feature might also relieve the need to 
place the ion producing foil in a reentrant cone, yielding a 



simpler target design ( Fernandez et al. 2009 1 . The esti 



mated parameters for a C ignitor beam are 400—500 MeV 



£ig includes only the energy of the proton beam. The total energy 
of the laser driver is Sig/rjp with rjp < 1 the conversion efficiency 
into protons. 
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energy per ion, > 10% efficiency, and 5£/£ < Q.2. To 



achieve such figures, mechanisms such as RPA (Sec IV. A) 



or BOA (Sec IV.C) might be more suitable than TNSA. 



Progress in related research has been recently reported 



by Hegelich et al. (2011) where separated experiments 
approaching the three above mentioned requirements are 
described. A fourth and so far unexplored issue might be 
the need to focus the ion beam. 

Another approach (Naumova et al. 2009) to ion- 



driven FI is based on RPA in the hole boring regime 



(Sec |IV.ArT ). In such scheme, differently from the above 
described ones, ion acceleration occurs in situ by direct 
interaction of an ultraintense, circularly polarized laser 
with the corona of the fusion plasma. |Tikhonchuk et af. 



(2010) report calculations on this scheme, assuming di- 
rect acceleration of Deuterons and characterizing possi- 
ble high gain regimes with ~ 12 — 17 kJ. This cor- 
responds to an overall ignition energy > 100 kJ and a 
required laser intensity exceeding 10^^ W cm~^. FI by 
laser-accelerated ions has been investigated theoretically 
in several other works|^ Integrated FI studies, on the 
route to ignition-class experiments, could be perfomed in 
either the electron or the ion approach in facilities to be 
developed, equipped with Petawatt-class laser systems. 



D. Biomedical applications 

Hadrontherapy is the radiotherapy technique that uses 
protons, neutrons or carbon ions to irradiate cancer tu- 
mors. The use of ion beams in cancer radiotherapjj^ 
exploits the advantageous energy deposition properties 
of ions as compared to more commonly used X-rays (see 
Figjs]): the range for a proton or ion is fixed by its energy, 
which avoids irradiation of healthy tissues at the rear side 
of the tumor, while the well localized Bragg peak leads to 
a substantial increase of the irradiation dose in the vicin- 
ity of the stopping point. The proton energy window of 
therapeutical interest ranges between 60 and 250 MeV, 
depending on the location of the tumor (the required 
Carbon ion range extends up to 400 MeV/nucleon). The 
typical dose of a treatment session is in the 1-5 Gray 
range, and typical currents are 10 nA for protons and 1.2 
nA for singly charged Carbon ions. 

Ion beam therapy has proven to be effective and ad- 
vantageous in a number of tumours and several clinical 



See e.g. Badziak et al.\ l|201lt; |Barriga-C arrasco et al.\ j2004|l 



IHoss cini Motlagh e t aL| | |2008[ l; Ramia and Ramrez| | |2004f ; |Shma^ 
I tov[ (^2003 2008 20TT]| 

" I'he use of energetic protons in radiology was first proposed by 
|Wilson| ([19461) and demonstrated by [Lawrence] | [1957[ l . Recent 
reviews on tfic state of tlie art in ion beam tfierapy are given, 
e.g., by [Amaldi and Kraft[ ( [2005[ l; [Scliardt et al.\ ([2010|l; and 
[Smithj | |2009[ |. Focus on researcli and possible improvements in 
therapy with heavy ions is given e.g. by [Kraft and Kraft| | j2009| l. 



facilities, employing mainly protons from synchrotron, 
cyclotron or linac accelerators are operational and rou- 
tinely treating a significant number of patients. While 
protons are the most widespread form of ion treatment, 
facilities using carbon ions also exist, as their higher bio- 
logical effectiveness makes them very effective in treating 



radioresistant and hypoxic tumours ( Schardt 2007 1 . In 
treatment centers, magnetic steering systems (Gantries) 
are employed for multi-directional irradiation of a laying 
patient. Gantries are costly, large and massive, with a 
weight exceeding 100 tons for proton systems, and 500 



tons for Carbon systems (Enghardt et al. 2011) 



The use of laser based accelerators has been proposed 
by several authors as an alternative to RF accelerators 
in proton and ion therapy system s ([Bulanov et al 



2002 



Bulanov and Khoroshkw| [2002{ [Fourkal et al.\ [2003b 



Malka et al. 2004), with potential advantages in terms 
of compactness and costs. Proposed options range from 
using laser-driven protons as high quality injectors in a 



RF accelerator ( Antici et al. 2011 1 to all-optical systems. 



in which the ion beam acceleration takes place in the 
treatment room itself and ion beam transport and deliv- 

It 



20021 



ery issues are thus minimized (Bulanov et al. 
is recognized that there are significant challenges ahead 
before laser-driven ion beams meet therapeutic specifi- 
cations, both in terms of maximum energy, energy spec- 
trum, repetition rate and general reliability, to the lev- 
els required by the medical and therapeutic standards. 



as critically reviewed by Linz and Alonso (2007) where 
specific issues are mentioned and a comparison with ex- 
isting accelerator technologies is made. At present, the 
ion beam parameters are still far from the requirements 
and it is clear that an extensive, long term activity will 
be needed to ascertain if and how laser-driven ion beams 
may become a competitive option. Several projects are 
currently active worldwide to explore the potential of 
laser-driven proton and ion sources for biomedical ap- 
plications, see e.g. Bolton et al. (2010); Borghesi et al. 



(2011); and Enghardt et al. (2011). In view of future ap- 



plications, several authors have started to design possible 
delivery systems, including target chamber and shielding 



( Ma et al. 2006 ) , particle energy selection and beam col- 



limation systems to enable operation with the broadband 



and diverging laser-driven beams (Fourkal et al. 


Hofmann et al. 


2011 


Nishiuchi et al. 


2010b). 



While currently a relative energy spread A£/£ ~ 10~^ 
is required for optimal dose delivery over the tumour re- 
gion, authors have also modelled approaches in which the 
native broad spectrum of laser-accelerated ions is used 
to obtain directly the Spread Out Bragg peak distribu- 
tions which are normally used to cover the tumour re- 



gion (Fourkal et al. 2007 Luo et al. 2008) and more 



in general advanced methods exploiting the properties of 



laser-accelerated beams (Schell and Wilkens 20101 



An important step in view of their future medical use 
is to assess the biological effect of laser-driven ions, and 
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highlight any anomaly associated to their pulsed 
trashort temporal profile 



ul- 



( Fritzler et al. 2003 Lefebvre et al. 2006 ) and econom 



Yogo et al. (2009) have first ically competitive with existing technologies. 



demonstrated the feasibility of cell irradiation studies 
using laser-driven protons, employing a suitable beam 
transport set-up and then applied a refined technique 
to infer, via a clonogenic assay, the Relative Biological 
Effectiveness (RBE) of ^ 2 MeV laser-accelerated pro- 
tons, as compared to irradiation with a standard X-ray 
source ( Yogo et al] 2011), in human cancer cells. The 
RBE observed (1.2 ± 0.1) is comparable with literature 
results employing RF-accelerated protons of comparable 
Linear Energy Transfer (LET) ( |Folkard[ [T996| . The dose 
required to cause significant cell damage (typically 1 to 
several Gray) was obtained in several irradiations tak- 



ing place at 1 Hz repetition rate. Kraft et al. (2010) also 



carried out proton irradiations of cells, highlighting dose- 
dependent incidence of double strand DNA break in the 



cells (Figj36). 

The peculiar characteristics of laser-driven protons 
have required the development of innovative dosimetric 
approaches, as described for example by [Fiorini et al.\ 
(2011) and Richter et al. (2011). In all these experi- 



ments, the dose (1-10 Gy) is delivered to the cells in 
short bursts of ~ ns duration. In experiments by [Kraft 



et al. (2010) and Yogo et al. (20111 the dose is fraction- 



ated and the average dose rate is comparable to the one 
used in irradiations with conventional acce lerator sources 
(~ 0.1 Gy s~^). In a recent experiment (Fiorini et al. 



2011) employing a high energy ps laser system, it has 



been possible to reach up to 5 Gy in a single exposure, 
reaching dose rates as high as 10^ Gy s~^. This allows 
access to a virtually unexplored regime of radiobiology. 



where, in principle, nonlinear collective effects (Fourkal 



et al. 2011 ) on the cell due to the high proton density in 



the bunch may become relevant. 

Besides cancer therapy, application of laser-driven ion 
beams in medical diagnosis has also been proposed. 
Multi-MeV proton beams can induce nuclear reactions 



E. Nuclear and particle physics 

The interaction of laser-driven high-energy ions with 
secondary targets can initiate nuclear reactions of vari- 
ous type, which, as mentioned before (Sec |ILE ) can been 
used as a tool to diagnose the beam properties. This also 
presents the opportunity of carrying out nuclear physics 
experiments in laser laboratories rather than in accelera- 
tor or reactor facilities, and to apply the products of the 
reaction processes in several areas. Reactions initiated 
by laser-accelerated high-Z ions have been studied in a 
number of experiments. McKenna et al. ( 2003a|b 2004) 
have shown that fusion reactions between fast heavy ions 
from a laser-produced plasma and stationary atoms in an 
adjacent "activation" sample create compound nuclei in 
excited states, which de-excite through the evaporation of 
protons, neutrons, and a-particles. A similar experiment 
with protons driving nuclear reactions and excitations in 



a Cu target has been reported by Hannachi et al. ( 2007 1 . 



Nuclear reactions of interest for spallation physics have 
also been investigated by employing the multi-MeV pro- 



ton beams (McKenna et al. 2005). The broad energy 
distribution of the beams is in this case advantageous for 
the determination of residual nuclide generation arising 
from specific spallation processes such as evaporation. In 
addition, MeV proton interaction with low-Z materials 
can produce short-lived isotope j^ o f medical interest, 
e.g. for PET diagnostic (see Sec V.D). 



Recently, a scheme of a "fission-fusion" process driven 



by RPA (Sec IV. A ) has also been proposed to produce 
neutron-rich nuclei in the range of the r-process (Habil 



et al. 2011); such studies, of relevance for astrophysi- 



in \ow-Z materials (SecV.E) in order to produce neu- 
trons, of possible interest for Boron Neutron Capture 
Therapy for cancer, or short-lived positron emitting iso- 
topes which may be employed in Positron Emission to- 
mography (PET) . PET has proven to be extremely useful 
in medical imaging of blood flow and amino acid trans- 
port and in the detection of tumors. Usually, reactions 
for PET are carried out by using up to 20 MeV protons 
or similar energy deuterons from cyclotrons with the con- 
comitant problems of large size and cost and extensive ra- 
diation shielding. Production of short-lived isotopes via 
laser-driven proton beams may be feasible in the near fu- 
ture with the possibility of employing moderate energy, 
ultrashort, high-repetition table-top lasers. Extrapola- 
tions based on present results point to the possibility of 
reaching the GBq activities required for PET if laser sys- 
tems capable of delivering 1 J, 30 fs pulses focused at 
lO^o W cm^^ with kHz repetition will become available L 



cal nucleosynthesis, would require intensities above 3 x 
10^2 w cm^^ that may be available with next-generation 
short-pulse laser facilities. 

Neutrons are an important product of the nuclear re- 
actions produced above, with potential applications in 
cancer therapy, neutron radiography, radiation damage 
of materials and transmutation of nuclear waste. The 
potential for laser-driven neutron sources is consider- 
able and offers advantages over accelerator- and reactor- 
driven sources in term of cost, compactness, brightness 
and short duration for applications such as fast neutron 



radiography ( Lancaster et al. . 2004 ) and studies of im 



pulsive damage of matter ( Perkins et al. 2000 ) . This has 



' Typical short lived positron emitters include ^^B, ^^C, ^-^N, ^^O, 
and ^®F. Related experiments have been reported e.g. by 
ICIarke et "aT] (|2006); Fritzler ef aT] ||2003|l; |Fujimoto ef aL[||2008| 
^9); Ledin gham et aL| | |2003| l; |Nemoto et a<.| ( |2001| l;"and Ogura| 



et al. 
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motivated several experiments on the production of neu- 
trons initiated by laser-driven proton beams on secondary 
targets. Experiments carried out at the VULCAN laser 



more general discussion of the required laser develop- 



facility have revealed neutron yields up to 4 x 10^ sr 
pulse at a laser intensity of 3 x 10^° W cm~^ (Yang et al. 



per 



2004a| ), produced via the "B(p,n)"C and 'U{ p,n)''Be 
reactions. The latter was also investigated by [Youssef 
et al. ( 2006 1 as a diagnostic of proton acceleration. Neu- 



tron production has also been observed in interaction 
with solid targets containing deuterium (typically deuter- 
ated plastic), which can be either directly irradiated by 



high intensity laser pulses (Disdier et al. 


1999 


Habara 


et al.\ 


2003 2004b||Norreys et al. 1998| Willinga 


e et al. 


2011a 


or irradiated by ions accelerated on a s 


eparate 
In both 


target 


(Fritzler et al. 2002 Karsch et al. 


2003). 



cases the neutrons are produced in the course of fusion re- 
actions of the type D(c?, n)'^He involving laser-accelerated 
deuterium ions as also observed in gaseous targets (Dit- 



mire et al. 1997 Grillon et al. 2002). Numerical model- 



ing and theoretical investigations of laser-driven neutron 
production have been carried out by several authors, see 



e.g. 


Davis and Petrov 


(2008 


2011 


); 


Ellison and Fuchs 


(2010 


); Macchi 


(2006 


); and 


Toupin et al. 


(2001 





The application of laser-accelerated ions in particle 
physics requires "by definition" the ions to be "relativis- 
tic", i.e. their total energy must exceed the rest en- 
ergy whose value per nucleon is ~ mp<? ~ 0.94 GeV. 
Presently, observed cut-off energies are more than an or- 
der of magnitude below this threshold. Nevertheless, the 
scalings inferred from either experiments or theoretical 
models and the foreseen availability of higher laser pow- 
ers in a few years suggests that GeV ions may eventually 
be produced and applied in selected particle physics ex- 
periments. Moreover, it may be noticed that the very 
low emittance that can be obtained for laser accelerated 
ions make them suitable for post-acceleration, e.g. as an 



injection source for heavy ion accelerators (Krushelnick 



et al. 2000j). Specific advantages might be the high num- 
ber of ions produced per shot combined with the short 
duration. 



Bychenkov et al. (2001b) estimated the threshold for 



production of pions by protons accelerated in a solid tar- 
get, obtaining that at intensities above 10^^ W cm~^ the 
flux of pions may be much higher than obtained with 
conventional accelerator techniques. It may be noticed 
that the prompt laser-driven, high field-gradient accel- 
eration of pions is of high interest because of the finite 
lifetime of such particles; a related discussion is reported 



b y|Mourou et a/.|(|2006 ). 

Pakhomov ( 2002 ) proposed the use of laser-accelerated 
protons at intensities of ~ 10^^ W cm~^ to drive, via pion 
generation and decay, pulses of 20 MeV muon-neutrinos 
of interest for, e.g., studies of neutrino oscillations. |Bu-| 
lanov et al. ( 2005 ) further explored this concept suggest- 



ing the radiation pressure dominated acceleration regime 



ments is reported by Terranova et al. ( 2006 1 



VI. CONCLUSIONS AND OUTLOOK 

We have reviewed about twelve years of research on 
ion acceleration driven by superintense laser pulses. The 
most investigated acceleration mechanism, namely the 
so-called Target Normal Sheath Acceleration, has been 
extensively discussed, surveying both the main experi- 
mental achievements and the underlying theoretical in- 
terpretation. In addition, we have provided an outlook 
to other proposed acceleration schemes, like Radiation 
Pressure Acceleration, Shock acceleration. Break-out af- 
terburner, acceleration in near-critical and underdense 
plasmas and resistively enhanced acceleration; for these 
mechanisms, the fundamental theoretical ideas have been 
presented, together with the most promising experimen- 
tal results. A brief presentation of the most promising 
possible applications of the laser-generated ion beams 
has been finally given. While we were completing our 
work, another extended review on this topic has been 



published (Daido et al. 2012), which we are delighted to 



(Sec IV. A. 2) as suitable for this class of experiments. A 



recommend as a complementary reading. 

This field of research has attracted an enormous inter- 
est and has shown unique potential for both innovative 
investigations and for applicative purposes. The impres- 
sive development in laser technology and the increasing 
use of advanced methods of material science for target 
manufacturing has resulted in a high level of sophisti- 
cation of current experiments, with new physical issues 
continuously emerging from the experimental investiga- 
tions. At the same time, theory and simulation naturally 
have found a fertile field, which poses original problems 
and suggest unexplored paths for reaching their solution. 
The vitality of this research area is well demonstrated 
by the need for frequent updates during the preparation 
of this Review, as new significant results continuously 
appeared in the literature. 

As discussed in the Introduction, future developments 
and achievements are naturally linked to foreseen devel- 
opments in laser technology, providing for the first time 
extreme laser intensities close to, or even beyond the limit 
at which the ions become relativistic. Emerging laser 
projects will enable to verify the scaling of ion acceler- 
ation physics in the ultra-high intensity regimes and its 
suitability for proposed applications, as well as to test 
theoretical ideas and provide to fundamental physics an 
example of "relativity in action" in a macroscopic, many- 
body system. The expected progress is not exclusively 
related to further developments of the solid state laser 
technology which has been so far the preferred route to 
producing high intensity laser pulses, as for example, re- 
cent experiments have attracted novel interest also in 
old technology such as C02 laser systems. Smart and 
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advanced target engineering, e.g. development of mul- 
tilayer, low-density, micro- and nanostructured targets, 
will also play an important role in establishing future 
directions of ion acceleration. 

The successful developments in this first period leads 
to a series of open questions, which will inform research 
in this field over the coming years. Will it be possible 
to reach and break the GeV/nucleon threshold? Will re- 
searchers achieve an active and satisfactory experimental 
control on the physics of laser-ion acceleration working 
on the detailed properties of laser pulses and of target 
material and structure? Is there hope to pursue front- 
edge nuclear and particle physics research in small scale 
laboratories thanks to the use of laser-driven ion beams? 
Will the research on ion acceleration result in practical, 
technological applications of direct societal benefit? We 
leave to the future experimental and theoretical research 
the answer to these and many other questions, some of 
which we probably do not even imagine today. 
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FIG. 25 a) Energy per nucleon vs. time from the analytical 
solution Eq.(32l of the LS model with R = 1. The dashed 
line gives the asymptotic ~ t^^'^ behavior, b) Scaling of the 
energy per nucleon as a function of the dimensionless pulse 
fluence o^t (where r is the pulse duration in units of the laser 
period) and of the surface density C [Eq.Q] for C = 1 (black 
line), 3.16 (green), 10 (blue), 31.6 (orange) and 100 (red). 
The values on the upper horizontal axis give the fluence in 
J cm~^ corresponding to Oqt for A — 0.8 /im. 
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FIG. 26 Three-dimensional simulations of thin foil accelera- 
tion in the Radiation Pressure Dominant regime (see text for 
parameters). Top frame: snapshots at f = 40r of ion density 
isosurface and Poynting vector in the y = plane. Bottom 
frame: The maximum ion kinetic energy versus time and the 
ion phase space projection (a;,Pa;) at f = SOT. The solid line 
corresponds to the analytical calculation according to the LS 
model. Reprinted figure with permission from'Esirkep ov et al\ 
12004), Phys. Rev. Lett. 92, 175003. Copyright 2013 by the 
American Physical Society. 
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FIG. 27 a) Three typical narrow band spectra for ZjA = 1/2 
impurity ions observed from thin (0.1 ptm) metallic targets 
(for parameters of shots 1-3 see explanation in part (b)). b) 
Peak energies for both ions with ZjA = 1/2 (filled squares) 
and protons with ZjA = 1 (empty squares) for seven shots 
with different polarization parameters [see e.g. Eq.(lll]; 
e = (LP), 0.47 (EP) and 0.88 (CP). The peak energy is 
plotted as a function of the parameter a^Tp/x, which corre- 
sponds to {mp/2me){yiTp) in our notations. The parameter 
set [ao, target material, thickness (/im), polarization] for the 
data points 17 is [15.5, Cu, 0.1, LP], [10, Cu, 0.05, CP], [13.8, 
Cu, 0.1, EP], [7.5, Al, 0.1, LP], [6.9, Al, 0.1, EP], [13.6, Al, 
0.5, CP], and [14.1, Al, 0.8, LP], respec tively. The red circl e 
marked with "H" is the data point from|Henig "eTari pOOQcl ). 
The solid line is the LS scaling (|35[). Reprinted figure with 

f2012") , Phys. Rev. Lett. 109, 
;he American Physical Society. 
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FIG. 28 Proton spectra from CO2 laser interaction with a 
Hydrogen gas jet ( [Haberberger et aL[|2012| ). Frame a) shows 
diff'erent spectra for a smooth (long) pulse (lower line) and a 
pulse train of 3 ps spikes (upper line); in the latter case, a peak 
appears in the spectrum. Frame b) shows narrow spectra ob- 
tained in different shots. See text for parameters. Reprinted 



by permission from Macmillan Publishers Ltd: (Haberberger 
et al 



20121, Nature Physics 8, 95. Copyright 2012. 
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FIG. 29 (Color online) Spectra of C®+ ions from laser in- 
teraction with ultrathin targets in the regime of relativistic 
transparency, as a function of target thickness and laser po- 
larization. Reprinted figure with permission from |Henig et aL\ 
(2009b|, Phys. Rev. Lett. 102, 095002. Copyright 2013 by 



the American Physical Society. 




FIG. 30 (Color online) Sketch of magnetic-field sustained ac- 
celeration of ions, showing the topology of the magnetic and 
electric fields and the flows of escaping and returning electron 
currents. Reprinted figure with permission from |Bulanov an"dl 
|Esirkepov| 12007] ), Phys. Rev. Lett. 98,049503. Copyright 
2013 by the American Physical Society. 
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FIG. 31 Proton probing of magnetic fields ( Cecchetti et al. 



20091. a), c): probing deflectogram of a laser- irradiated foil 
IS^ulse, lO^^^ W cm ^ on a 6 /im AI foil) obtained with a 
5.5 MeV proton beam facing the foil, entering from the non- 
irradiated side in a) and from the opposite, laser-irradiated 
side in c). The inversion of the deflection pattern reveals 
the effect of a toroidal B-field (the asymmetrical pattern in 
c) is due to a nonideal intensity distribution in the focus). 

b) , d): particle tracing simulations for the conditions of a), 

c) assuming a suitably parametrized B-field. Reprinted with 



permission from Physics of Plasmas: Cecchetti et al. (20091, 
Phys. Plasmas 16, 043102. Copyright 2009, American Insti- 
tute of Physics. 
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FIG. 32 (Color online) Left frame shows the proton imaging 
set-up with increased dynamical range in the time domain 
(Quinn et al. 2009a|b I. By placing a wire target at an an- 
gle with respect to the probe beam it is possible to resolve 
the propagation at a velocity close to c of a field front along 
the laser-irradiated wire, as shown in the right frames (top: 
experimental images, bottom: particle tracing simulations) . 
Reprinted figure with permission from Quinn et al.\ { 2009a I , 
Phys. Rev. Lett. 102, 194801. Copyright 2013 by the Amer- 
ican Physical Society. 
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FIG. 33 (Color online) Heating of solid targets by protons, 
(a) Experimental setup for flat and focusing target geome- 
tries. Each target consists of a flat or hemispherical fO fim 
thick Al foil irradiated by the laser, and a flat 10 thick 
Al foil to be heated by the protons, (b) Corresponding 
streak camera images showing space- and time-resolved ther- 
mal emission at 570 nm from the rear side of the proton- 
heated foil. Proton focusing by the hemispherical foil leads 
to a stronger, more localized heating. Reprinted figure with 
permission from Patel et al.\ ( 2003 1, Phys. Rev. Lett. 91, 
125004. Copyright 2013 by the American Physical Society. 




FIG. 34 (Color online) Set-up for X-ray probing of a solid 
Al target heated by protons: (a) top view of the experiment, 
(b) schematic of the x-ray spectrometer; (c, d) snapshots of 
temperature and density profiles of the heated 1.6 ^m thick 
Al foil as given by self-consistent simulations, demonstrating 
isochoric heating up to 1 5 eV temperatures. R eprinted fig- 
ure with permission from Mancic et aZ.| (2010al, Phys. Rev. 
Lett. 104, 035002. Copyright 2013 by the American Physical 
Society. 
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FIG. 35 Concept of proton-driven Fast Ignition in the TNSA- 
based, cone guided scheme (Key 20071. Typical parameters 



required for the ion beam and optimization issues are also 
indicated. Reprinted with permission from Physics of Plas- 
Key| ( |2007| ), Phys. Plasmas 14, 055502. Copyright 2007, 



American Institute of Physics. 
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FIG. 36 (Color online) Top: overview of an experimental 
setup for integrated dosimetry and cell irradiation system by 
laser accelerated protons. Bottom: fluorescence microscopy 
view of SKX tumour cell nuclei irradiated with such sys- 
tem, showing that the number of DNA double-strand breaks 
(bright spots, yellow-pink in the color version) increases with 
the delivered dose. Reprinted figure from Kraft et al. (20101, 
New J. Phys. 12, 085003. By permission from Institute of 
Physics Pubhshing (2013). 



